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ABSTRACT 

We study the star formation efficiency (SEE) in simulations and observations of turbulent, mag- 
netized, molecular clouds. We find that the probability density functions (PDFs) of the density and 
the column density in our simulations with solenoidal, mixed, and compressive forcing of the turbu- 
lence, sonic Mach numbers of 3-50, and magnetic fields in the super- to the trans- Alfvenic regime, all 
develop power-law tails of flattening slope with increasing SFE. The high-density tails of the PDFs 
are consistent with equivalent radial density profiles, p oc r~'^ with ~ 1.5-2.5, in agreement with 
observations. Studying velocity-size scalings, we find that all the simulations are consistent with 
the observed v ex scaling of supersonic turbulence, and seem to approach Kolmogorov turbulence 
with V oc below the sonic scale. The velocity-size scaling is, however, largely independent of 
the SFE. In contrast, the density-size and column density-size scalings are highly sensitive to star 
formation. We find that the power-law slope a of the density power spectrum, PsBip^k) cx /c", or 
equivalently the A- variance spectrum of the column density, (7^(11, cx switches sign from a <0 
for SFE ~ to a > when star formation proceeds (SFE > 0). We provide a relation to compute 
the SFE from a measurement of a. Studying the literature, we find values ranging from a = —1.6 to 
-hi. 6 in observations covering scales from the large-scale atomic medium, over cold molecular clouds, 
down to dense star-forming cores. From those a values, we infer SFEs and find good agreement with 
independent measurements based on young stellar object (YSO) counts, where available. Our SFE-a 
relation provides an independent estimate of the SFE based on the column density map of a cloud 
alone, without requiring a priori knowledge of star-formation activity or YSO counts. 
Subject headings: ISM: clouds - ISM: kinematics and dynamics - ISM: structure - magnetohydrody- 
namics (MHD) - stars: formation - turbulence 



1. INTRODUCTION 

The most important physical processes determining 
star formation are turbulenc e, gravity, a nd magnetic 
fields (see the reviews by Mac Low fc Klessen 2004; 



Elmegreen fc Scalo 2004) |Scalo k ]jlmegreenP004| |Mc- 
Kee &: Ostrikern2007j ). On one hand, molecular cloud 
turbulence — because it is supersonic — compresses inter- 
stellar gas in shocks and filaments. If a critical amount 
of mass is swept up in a local compression, that gas 
can become gravitationally unstable and collapse, giving 
birth to new stars. On the other hand, turbulence and 
magnetic fields together carry about the same amount 
of energy as there is in gravitational binding energy in a 
whole molecular cloud (e.g., Stahler & Palla 2004). Thus, 
both turbulence and magnetic fields provide large-scale 
support against global cloud collapse. That support is 
needed, because otherwise molecular clouds w ould simply 
collapse as a whole, which is not observed (Zuckerman 
Palmer II 19 74 Zuckerman & Evans 1974). I'he role of 



gravity and magnetic pressure is unidirectional: gravity 
always acts to promote collapse, while magnetic pressure 
always counteracts collapse. The role of turbulence on 
the other hand is dual: it provides global support, but 
at the same time produces local compression — the seeds 
for stellar birth. 
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Hence, turbulence, gravity, and magnetic fields are cru- 
cial for our understanding of star formation. These three 
physical effects are likely the major players in controlling 
the star formation rate (SFR) in the Milky Way and po- 
tentially in oth er galaxies, as discussed in jPeder rath fc| 



|Klessen| ( [2012 



hereafter Paper I). There, we compared 
theoreTical and numerical models of the SFR with obser- 
vations of Galactic clouds, finding very good agreement 
for a wide range of physical parameters, including ex- 
treme cases of turbulent forcing (solenoidal versus com- 
pressive forcing of the turbulence), sonic Mach numbers 
between 3 and 50, and magnetic fields ranging from the 
super- Alfvenic to the trans- Alfvenic regime, all in the 
range of parameters observed in real molecular clouds. 

Despite this large variety of physical parameters, the 
conversion of gas into stars is always quite inefficient in 
molecular clouds, with SFEs of typicall y on l y a few per 



cent (Myers et al.||1986[ Evans et al. 2009} |Lada et ah 
'2010). However, it is still poorly understood why this 



is the case. It is the aim of this study to shed light on 
this subject by providing a detailed analysis of the SFE 
in the simulations of Paper I, and comparing them to 
observations of Galactic clouds. 

The SFE of clouds is often uncertain or unknown. The 
main difficulty in estimating the SFE lies in obtaining a 
clear census of YSOs, in particular in dense and confused 
regions where many protostellar disks are embedded and 
closely packed, such that the YSO population is often 
underestimated in observations. There is some indica- 
tion that the SFE ranges from less than 1% up to 50% 
in most clouds, depending on the particular cloud stud- 
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ied, its physical parameters, the actual subregion con- 
sidered inside that cloud, and the evolutionary stage of 
that region. This is a very wide range of SFEs. The 
aim of this paper is to advance our understanding of the 
physical processes that determine this range and provide 
measures to constrain the SFE for a particular cloud re- 
gion observed. For instance, considered as a whole, giant 
molecular clouds typically only con vert at most a few 
percent of their g as into stars (e.g., Myers et al. 1986; 
[Evans et al.||2QQ9 ). In contrast, when looking at regions 
of dense cor e and cl uster formation, observers find higher 
Wilking fc Lada ||1983l IBeuther et al. 



SFEs 



^ ||1983| |Beu 

Koenig et al.||2Ud5 ), eventually approaching the local ef- 



2002; 



ticiency e on scales of individual protostellar cores. This 
local efficiency e determines the gas fraction of a single 
protostellar accretion envelope that actually ends up on 
the protostar(s) forming in a dense core. Although gas 
is efficiently falling onto the protostar, jets, winds and 
outffows launched from the protostellar disk counteract 
the inflow, thereby limiting accretion such that e < 1. 
We estimated this local efliciency with e ~ 30%-70% by 
comparing our simulatio ns with Galactic observations by 
Heiderman et al.| (|2010|) in Paper I. It is thus plausible 



that the SFE naturally approaches the local efliciency e 
on sufliciently small scales, when a single, dense, star- 
forming core is considered. 

Here, we investigate the dependence of the SFE on 
molecular cloud scales and on physical parameters, such 
as the driving of turbulence, the sonic and Alfven Mach 
numbers, and the virial parameter of a cloud. First, we 
study the structure of the column density and how it 
changes with SFE in Section^ We then discuss the de- 
pendences of the probability density functions (PDFs) of 
the volumetric and column densities in Section lU We 
find that the PDFs depend on the forcing of the tur- 
bulence and on the sonic and Alfven Mach numbers. 
The PDFs develop power-law tails of fiatting slope in 
all our numerical models when s tar formation proceeds 
as seen in observations of clouds (Kainulainen et al. 12009 



[Schneider et al. 2012), suggesting that the PDI^' can be 
used to infer the SFE. In Section [5) however, we find 
that a more reliable and independent distinction between 
star-forming and quiescent clouds is possible by study- 
ing the spatial scaling of column density with Fourier 
or A-variance analyses ( Ossenkopf et al.| 2001), i.e., by 
studying the density-size scaling in molecular clouds. We 
extend this qualitative discriminator of star formation 
to a quantitative method for estimating the SFE in Sec- 
tion [6| by measuring how the slope a of the density power 
spectrum changes with SFE, depending on the forcing of 
turbulence, the sonic Mach number, and the magnetic 
field. Inverting this relation yields a function SFE(a), 
which can be used to estimate SFE from a measurement 
of a in a dust or integrated molecular line map. The ad- 
vantage of this new method is that no a priori knowledge 
of star formation activity or YSO counts is required to 
estimate SFE. In Section [7| we apply this method to ob- 
servations and compare inferred SFEs with independent 
estimates, indicating good agreement. In Section we 
discuss uncertainties and limitations of the new method. 
Our conclusions are summarized in Section [9l 



2. NUMERICAL SIMULATIONS 



Our numerical simulation techniques are explained in 
detail in Paper I. Here we only give a brief overview of the 
most important aspects of the simulation methodology 
and provide a complete list of simulation parameters as 
in Paper I. 
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We use the ad aptive mesh refinement (AMR, Berger 
~ ' " T9891) code F LASH v2.5 dFryxeh et"^ 

model isothermal 



^ _^ 2008[ ) to 

gravitating, magnetohydrodynamic (MHD) turbulence 
on three-dimensional (3D), periodic grids with resolu- 
tions of TV^es = 128^-1024^ grid points. These are ah 
uniform- grid simulations, except for one N^^^ = 1024^ 
simulation, where we use a root grid with 512^ cells 
and one level of AMR with a refinement criterion to 
ensure that the local Jeans lengths is covered with at 
least 32 grid cells, in order to resolve turbulent vort icity 



and ma gnetic- field amplification on the Jeans scale (Sur 
et al.|2010^, Federrath et al.|2011c||Turk et al.|2012D . Tor 



solving the MHD equati ons, we usethe HLL 3K positive- 
definite Riemann solver (Waagan et al.|20"TT ). The MHD 
equations are closed with an isothermal equation of state, 
which is a reasonable approximation for dense, molecu- 
lar gas of solar metallicity, over a w ide range of densities 
(Wolfi re et al.||1995l [Omukai et al.|[20 05; Pavlovski et al. 
| 2006i [Glover ^ Mac Low|2007a|bl|GTQver et al |2010( |TTiTT 
|et al.poill [Hennemann et al.||2012p . The self-gravityot 
the gas is computed with a multi-grid Poisson solver (see 
|Ricker„2008| ). 

2.1. Turbulent Forcing 

To drive turbulence in the simulations, we apply a 
stochastic acceleration field Fgtir as a momentum and en- 
ergy source term. Fgtir only contains large-scale modes, 
1 < /c < 3, where most of the power is injected at the 
k = 2 mode in Fourier space, which corresponds to half 
of the box size L in physical space. We thus model tur- 
bulent forcing on large s cales, as favored by molecular 



cloud observations (e.g., Ossenkopf fc Mac Low |2002 
Heyer et al. 2006; B runt et a l. 2009 ; jKoman-Duval et af. 
2011D. Smaller scales, > 3, are not affected directly 



by tne forcing, such that turbulence can develop self- 
consistently there. We use the Ornstein-Uhlenbeck (OU) 
process to model Fgtir, which is a well-defined stochastic 
process w ith a finite autocorrelation timescale ( Eswaran] 
fc Pope| [19881 [Schmidt et"ar]|2006D . We set the au- 
tocorrelation time equal to the turbulent crossing time 
on the largest scales of the system, T = L/(2A^Cg), 
with the RMS sonic Mach number Ai and the constant 
sound speed Cg = 0.2kms~^, leading to a smoothly vary- 
ing stochastic forcing in space and time (for details, see 
Schmidt et al. 2009; Federrath et al. 2010b; Konstandm] 
et al. 2012a, Paper I). 

We can adjust the mixture of solenoidal and compres- 
sive modes of our turbulent forcing, by applying a pro- 
jection in Fourier space. Here we compare simulations 
with solenoidal forcing (V • Fgtir = 0) and compressive 
forcing (V x Fgtir = 0), as well as an intermediate mix- 
ture of both. Solenoidal and compressive forcing are ex- 
treme cases, while in real molecular clouds we expect 
some mixture with a range of possible ratios between 
solenoidal and compressive modes (see the discussion of 
physical drivers of turbulence and their m ode charac- 
teristics in Paper I). For instance, ,Motte et al. (1998) 
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describe a scenario in which the observed ahgnment of 
young protostars in the p Ophiuchi central region could 
have been triggered by a shock compression, induced by a 
supernova/ wind shell and/or expanding H ll regions from 
nearby OB associations. Such expanding shells are com- 
pressive forcing mechanisms, because they primarily ex- 
cite compressible velocity modes. 

2.2. Sink Particles 

In order to model collapse and accretion, we use an ad- 
vanced AMR-based approach for sink particles, in which 
only bound and collapsing gas is acc reted (for a detailed 
analysis, see ^ Federrath et aL||2QlQa ). The key feature of 
this approach is to detine a control volume centered on 



grid cells exceeding a density threshold. [Truelove et al. 
(l997) found that the Jeans length must be resolved with 
at least 4 grid cells to avoid artificial fragmentation, lead- 
ing to a resolution-dependent density threshold criterion 
for triggering sink particle creation checks. 



Psink 



(1) 



where the sink particle accretion radius rgink is set to 
2.5 grid-cell lengths at the maximum level of refinement, 
corresponding to half a Jeans length at Psink- This guar- 
antees that the Jeans length is still resolved with 5 grid 
cells prior to potential sink particle creation to avoid ar- 
tificial fragmentation. 

Grid cells exceeding the density threshold given by 
Equation ([T]), however, do not form sink particles right 
away. First, a spherical control volume with radius rgink 
is defined around the cell exceeding Psink, in which a se- 
ries of checks for gravitational instability and collapse are 
performed. If all checks are passed, a sink particle is cre- 
ated in the center of the control volume. This procedure 
avoids spurious sink particle formation, and allows us to 
trace only truly collapsing and star-forming gas. 

Once a sink particle is created, it can gain mass by 
accreting gas from the AMR grid, but only if this gas 
exceeds the threshold density, is inside the sink particle 
accretion radius, is bound to the particle, and is collaps- 
ing toward it. If all these criteria are fulfilled, the excess 
mass above the density threshold defined by Equation ([l]) 
is removed from the MHD system and added to the sink 
particle, such that mass, momentum and a ngular mo- 
mentum are conse rved by construction (see Federrath] 



et al. 



2010a', 'MTa, for details). 



Sink particle gravity is computed by direct A^-body 
summation over all sink particles and grid cells. We use 
a second-order accurate Leapfrog integrator to advance 
the sink particles on a timestep that allows us to resolve 
close and highly eccentric orbits of sink particles with- 
out introducing sign ificant errors on super-r esolution grid 
scales, as tested in [Federrath et al. (2010a). 



2.3. Initial Conditions, Procedures, and List of Models 

We start our numerical experiments with uniform den- 
sity and zero velocities. The forcing term Fgtir (see Sec- 
tion |2.1[ ) drives random motions, until a state of fully- 
developed, supersonic turbulence is reached after two 
large-scale turbulent crossing tim es, 2T = LHJyic^), 
as found in previous studies (e.g . , [Kl essen et al." 2000; 
[Klessen 2001, Heitsch et al.|[2QQTl 



2010b| [Price fc Federrat"h||2010[ picic et al.|[2QT2 |. This 
lully-aeveloped turbulent state is the initial condition 
for our star-formation experiments, when self-gravity is 
added and formation of sink particles is allowed. We 
denote this time as t = in the following. 

All our numerical simulations and their basic param- 
eters are listed in Table [l] These parameters were cho- 
sen to roughly follow observed properties of molecular 
clouds, covering a range of cloud sizes L ~ 0.3-200 pc, 
masses Mc ~ 300 to 4 x 10^ Mq, and velocit y disper- 
sions cFy ^ 0.6-10 km s~^ (e.g., Larson 198l| [Solomon 
et al 119871 IFalgarone et al. 1992; Ossenkopf & Mac Low 
^QOWlHeyer Brunt||2004, .Heyer et al...2009, .Roman- 



et al 119871 IFalgarone et al. 1992; Ossenkopf & Mac Low 
20021 IHeyer fc Brunt||2004, .Heyer et al..,2009, .Roman- 
Duval et al. 201 ip , wi th typical cloud scalings sum ma- 
rized and discussed in M ac Low fc Kle ssen[ ( [2004 ) and 



[McKee fc Ostriker| ([2007|). For the MHD models, we 
chose magnetic field strengths consistent with the r ange 
observed in clouds of that size and density (see Crutcher 
19991 [Heiles fc Troland[[2QQ5|)Crutcher et al.[|2()10| ). We 



compare models with initial line-ot-sight magnetic field 
strengths Bq^z = 1, 3, and 10 /iG. 

2.4. Definition of the SFE 



We use t he standard definition of the SFE (e.g., Myers 
et al.|1 986^), which is the mass in star- forming gas (stars), 
M^, divided by the total mass of the cloud. 



SEE 



Mr' 



(2) 



We measure M^ in the simulations simply by taking the 
sum of all sink particle masses formed in the simulation 
at a given time. The denominator in Equation ([2| is 
constant over time and given by the total cloud mass, 
Mc, in Table [T] because our computational boxes are 
closed, without global gas infiow or outfiow. Gas that 
is turned into sinks is removed from the gas phase to 
conserve the total mass in the system. 

3. EVOLUTION OF THE CLOUD STRUCTURE 
After the initial turbulent state has been established 



(see Section 2.3), we study the subsequent evolution un- 



Terrath et al. 2009, 



der the infiuence of self-gravity. First, we look at how the 
column density structure changes with increasing SFE 
in Figures [T] and [2] We see that density contrasts in- 
crease with increasing sonic Mach number (from left to 
right: M ^ 10, and 50). For a given Mach number, 
solenoidal forcing (Figure IT]) produces significantly less 
variation in the column density than compressive forc- 
ing (Figure J2|. Star formation (increasing SEE from top 
to bottom m the figures) affects mostly the small scale 
structure where collapse proceeds, but also reaches to 
large scales, best seen in the runs with ^ 3. There, 
we also find evidence for global collapse, while star for- 
mation is locally much faster in the high-Mach number 
{M. ~ 10 and 50) cases, such that no significant global 
collapse occurs even when SFE = 20% is reached. This 
is because in the high-Mach number cases, turbulence 
alone provides significant local compression, such that 
small-scale structures quickly proceed to collapse. For 
instance, in the most extreme case {M. ~ 50 with com- 
pressive forcing; right panels in Figure [2|, the large-scale 
structure does not change at all during the short time 
required to wind up 20% of the gas in bound cores and 
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Table 1 

Basic Parameters of the Numerical Models of Forced, Supersonic, Self-gravitating, (Magneto)hydrodynamic Turbulence. 



Model 




A^res 


Forcing 


po [g cm- 




L[pc] 


Me [Mo] 


ay [km/s] 










M 


b 




Ma 


(1) 




(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


01) GT256sM3 




256 


sol 


5.8x10- 


-19 


3.3x10-1 


3.1x10^ 


0.59 





oo 


0.07 


1.4 


2.9 


1/3 


oo 


oo 


02) GT512sM3 




512 


sol 


5.8x10- 


-19 


3.3x10-1 


3.1 xlO^ 


0.59 





oo 


0.07 


1.4 


3.0 


1/3 


oo 


oo 


03) GT256mM3 




256 


mix 


5.8x10- 


-19 


3.3x10-1 


3.1 xlO^ 


0.61 





oo 


0.08 


1.1 


3.1 


0.4 


oo 


oo 


04) GT256cM3 




256 


comp 


5.8x10- 


-19 


3.3x10-1 


3.1x10^ 


0.58 





oo 


0.07 


0.46 


2.9 


1 


oo 


oo 


05) GT512cM3 




512 


comp 


5.8x10- 


-19 


3.3x10-1 


3.1 xlO^ 


0.58 





oo 


0.07 


0.48 


2.9 


1 


oo 


oo 


06) GT256sM5 




256 


sol 


3.3x10- 


-21 


2.0x10° 


3.9 xlO^ 


1.0 





oo 


1.0 


8.0 


5.0 


1/3 


oo 


oo 


07) GT256mM5 




256 


mix 


3.3x10- 


-21 


2.0x10° 


3.9 xlO^ 


0.99 





oo 


0.98 


5.4 


5.0 


0.4 


oo 


oo 


08) GT256cM5 




256 


comp 


3.3x10- 


-21 


2.0x10° 


3.9x10^ 


0.91 





oo 


0.82 


1.5 


4.5 


1 


oo 


oo 


09) GT128sM10 




128 


sol 


8.2x10- 


-22 


8.0x10° 


6.2x10^ 


2.1 





oo 


1.1 


11. 


10. 


1/3 


oo 


oo 


10) GT256sM10 




256 


sol 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


2.1 





oo 


1.1 


12. 


10. 


1/3 


oo 


oo 


11) GT512sM10 




512 


sol 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


2.1 





oo 


1.1 


12. 


10. 


1/3 


oo 


oo 


12) GT512mM10 


(si) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2x10^ 


2.1 





oo 


1.1 


4.5 


11. 


0.4 


oo 


oo 


13) GT512mM10Bl 


(si) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


2.1 


1 


8.2 


1.1 


5.4 


10. 


0.4 


2.8 


12. 


14) GT512mM10 


(s2) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


2.2 





oo 


1.2 


8.4 


11. 


0.4 


oo 


oo 


15) GT512mM10Bl 


(s2) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2x10^ 


2.2 


1 


8.2 


1.2 


9.5 


11. 


0.4 


1.8 


10. 


16) GT256mM10 


(s3) 


256 


mix 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


2.0 





oo 


1.0 


5.9 


10. 


0.4 


oo 


oo 


17) GT512mM10 


(s3) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


2.0 





oo 


1.0 


5.9 


10. 


0.4 


oo 


oo 


18) GT512mM10Bl 


(s3) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2x10^ 


2.0 


1 


8.2 


0.97 


6.4 


9.9 


0.4 


3.6 


13. 


19) GT256mM10B3 


(s3) 


256 


mix 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


1.8 


3 


0.92 


0.81 


8.4 


9.0 


0.4 


0.20 


2.9 


20) GT512mM10B3 


(s3) 


512 


mix 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


1.8 


3 


0.92 


0.83 


8.7 


9.1 


0.4 


0.18 


2.7 


21) GT256mM10B10 


(s3) 


256 


mix 


8.2x10- 


-22 


8.0x10° 


6.2x10^ 


1.8 


10 


0.08 


0.79 


6.6 


8.9 


0.4 


0.04 


1.3 


22) GT128cM10 




128 


comp 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


1.8 





oo 


0.81 


1.2 


9.0 


1 


oo 


oo 


23) GT256cM10 




256 


comp 


8.2x10- 


-22 


8.0x10° 


6.2 xlO^ 


1.8 





oo 


0.85 


1.1 


9.2 


1 


oo 


oo 


24) GT512cM10 




512 


comp 


8.2x10- 


-22 


8.0x10° 


6.2x10^ 


1.9 





oo 


0.87 


1.1 


9.4 


1 


oo 


oo 


25) GT256sM20 




256 


sol 


2.1x10- 


-22 


3.2x101 


9.9 xlO^ 


4.1 





oo 


1.0 


11. 


20. 


1/3 


oo 


oo 


26) GT256mM20 




256 


mix 


2.1x10- 


-22 


3.2xl0i 


9.9x10^ 


4.2 





oo 


1.1 


4.5 


21. 


0.4 


oo 


oo 


27) GT256cM20 




256 


comp 


2.1x10- 


-22 


3.2xl0i 


9.9 xlO^ 


4.0 





oo 


1.0 


0.60 


20. 


1 


oo 


oo 


28) GT256sM50 




256 


sol 


3.3x10- 


-23 


2.0x10^ 


3.9x10^ 


10. 





oo 


1.1 


12. 


52. 


1/3 


oo 


oo 


29) GT512sM50 




512 


sol 


3.3x10- 


-23 


2.0x10^ 


3.9 xlO^ 


10. 





oo 


1.1 


13. 


52. 


1/3 


oo 


oo 


30) GT256mM50 




256 


mix 


3.3x10- 


-23 


2.0x10^ 


3.9x10^ 


10. 





oo 


1.0 


7.0 


51. 


0.4 


oo 


oo 


31) GT512mM50 




512 


mix 


3.3x10- 


-23 


2.0x10^ 


3.9 xlO^ 


10. 





oo 


1.1 


7.4 


51. 


0.4 


oo 


oo 


32) GT256cM50 




256 


comp 


3.3x10- 


-23 


2.0x10^ 


3.9 xlO^ 


9.8 





oo 


0.95 


0.54 


49. 


1 


oo 


oo 


33) GT512cM50 




512 


comp 


3.3x10- 


-23 


2.0 xlO^ 


3.9 xlO^ 


9.9 





oo 


0.99 


0.56 


50. 


1 


oo 


oo 


34) GT1024cM50 




1024 


comp 


3.3x10- 


-23 


2.0x10^ 


3.9x10^ 


10. 





oo 


1.00 


0.55 


50. 


1 


oo 


oo 



Notes. Column (1): simulation name. Columns (2-10): maximum grid resolution A^res in one direction of the 3D, cubic domain, mode 
of forcing (solenoidal, mixed, compressive), mean density po, linear box size L, total mass Mc, velocity dispersion on the box scale cry, 
mean magnetic field strength Bq (in z-direction of the domain), initial plasma /5o, and virial parameter Q;vir,o = 5a^L/(6GMc) for a 
spherical-cloud approximation with uniform density. Columns (11-15): time-averaged virial parameter oivir = 2£'kin/ l^'gravl of the 3D, 
inhomogeneous density field in the regime of fully-developed turbulence, sonic Mach number forcing parameter 6, ratio of thermal to 
magnetic pressure (plasma (3), and Alfven Mach number A4a- To guide the eye, horizontal lines separate models with different sonic Mach 
number. See Paper I for details of the simulation parameters. 

stars. In contrast, the turbulent density seeds in the 
J\4 ^ 3 runs have such a low amplitude that large-scale 
gravitational contraction is required to increase the local 
densities to the point when star formation can proceed. 
This behavior is reflected in the SFRs of all our simula- 
tions, ranging over two orders of magnitude, depending 
on the virial parameter, the turbulent forcing, the sonic 
Mach number, and the Alfven Mach number, studied in 
detail in Paper I. 

The column density evolution and dependence on the 
forcing and the sonic Mach number give us a visual im- 
pression of the time dependence and spatial distribution 
of star formation. We immediately see that the density 
variance depends on the sonic Mach number and the forc- 
ing already in the purely turbulent regime (top rows in 
Figures [l] and [2| and become stronger on small scales, 
where star formation occurs. Stars form primarily in 
dense filaments, with more massive stars being located 



mostly at intersections of multiple filaments or where fil- 
aments seem to bend. There, mass is transported to the 
protostars more efficiently, because it is funneled through 
multiple filaments toward their intersections where grav - 
itational focusing occurs (Burkert & Hartmann 2004). 
Similar structures and the correlation of star formation 
with filaments and their in tersections have been observed 
in molecular clouds (e.g., Andre et al.| 2010; Schne ider 
et al.|[2QTQ{ [Arzoumanian et al. | |:^011r[Schneider et a 



2012). Supersonic turbulence and gravity are the nat- 



ural physical processes that can produce such filaments 
in gas that can cool sufficiently to roughly maintain a 
constant temperature over a wide range of densities (see 
^Peters et al. 20 12, and references therein), as seen in the 
observations and in the present simulations. 

In order to gain more insight into the correlation of gas 
density with star formation, we investigate density and 
column density PDFs in the following. 
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Moch 3, sol, Bq = 
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il/box = 3.9E+06 Mq 
^srnk= L,,, = 200.0 pc 
SFE= 0.0% r., = 



SFE= 0.0% 



t/t„ = 0.66 Moch 50, sol, Bq = t/t„ = 0.25 





JI/,,, = 3.lE+02 Mo 
^box= 0.3 pc 
r3;„, = 335 AU 
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i M^, = 6.2E+03 Mo 
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-2,6 



-2,8 



-3,0 



Figure 1. Time evolution of column density (normalized to the mean column density Eq) for simulations with solenoidal forcing at 512^ 
resolution and ~ 3 (left), ~ 10 (middle), and ~ 50 (right), when t = (state of fully developed turbulence; top row), and when 
SFE = (right before the first sink particle forms; 2. row), SFE = 5% (3. row), and SFE = 20% (bottom row). Sink particles are shown 
as circles with radius rgink- Simulation parameters and variables are indicated in each panel and in Table ^ 
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^sln.= 91 

SFE = 20.3 % 




Figure 2. Same as Figure [l] but for compressive forcing of the turbulence. 
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4. THE DENSITY PDF 

The density PDF of the turbulent gas is a key for an- 
alytic model s of the core and stellar initi al mass func- 
tions (Padoan & Nordlund 2002; Hennebelle & Chabrier 
2008, 2009; Elmegreen 2011; Veltchev et al. 2011; Don 
et al. 2012; Parravano et al. 2012; Hopkins 2012c a); 



kov 
the 



Kennicutt- Schmidt relation "TKrumholz & McKee 2005 

TT, nrTTTTTzn ^^T^ rm ^ . . j . ■ < ^ 



|Tassis"2007|), the SFE 



and the^FT 



KrumhoIi"& McKee 2005 Padoan & JNordlund 2011 



Hennebelle & Chabrier 201 1| Paper I). We investigate 
the dependences of the PDF on the forcing of the tur- 
bulence, the Mach number, and the magnetic field, and 
pay specific attention to its dependence on the SFE in 
the following. 



4.1. The Density PDF of Supersonic, Isothermal 
Turbulence 
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Studies of non-gravitating, isothermal turbulence show 
that the density PDF is close to a log-normal distribu- 
tion, 

1 / {s-sor 



Pais) = 



exp 



expressed in terms of the logarithmic density, 
s = In (p/po) • 



(3) 



(4) 



The PDF is a normal (Gaussian) distribution in s, and 
hence a log-normal distribution in p. The quantities po 
and So denote the mean density and mean logarithmic 
density, the latter of which is related to the standard 
deviation as by sq = — due to the normalization 
and mass- conservation constraints of the PDF ( Vazquez- 
|Semadeni| |1994l |Li et al"][2QQ3g Federrath e t"l^ . 

The standard deviation (Js m Equation Q is a measure 
of how much the density varies in a turbulent medium 
and depends on 1) the amount of compression induced 
by the turbulent forcing mechanism, 2) t he sonic Mach 



number, and 3) the magnetic pressure. Molina et al. 
(2012) provide a rigorous derivation of the variance ot 
thiTT)F, 



cr? = In 1 + b^M^' 



(5) 



with the forcing parameter 6, the sonic RMS Mach num- 
ber A^, and the ratio of thermal to magnetic pressure, 
which can be expressed as the ratio of sound to Alfven 
speed or as the ratio of Alfven to sonic Mach number, 
P = 2cl/v\ = 2M\/M^. The forcing parameter h 
in Equation ([5| was shown to vary smoothly between 
h ~ 1/3 for purely solenoidal (divergence-free) forcing, 
and 6 ~ 1 for purely compressive (curly-free) for cing of 



the turbulence (see Section |2.1l Federrath et al. 2QQ8b 



Schmidt et al.|2QQ9l|FederratTret al .2QlQb,.Seitried et a. 



20111 [Micic et al.||2Q12| IKonstandin et al.|T2Q12ap . A 

stochastic mixture of fo rcm^ modes in 3D space leads to 
h - 0.4 (see Figure 8 in|Federrath et al.|[2010b). 

When gravity is included and significant collapse sets 
in, the d ensity PDF de velops a power-law tail at high 
densities ( |Klessen||2000 ), which we concentrate on in the 
following. 

4.2. The Density PDF of Self- gravitating Turbulence 
4.2.1. Volumetric Density PDFs 

Figure |3] shows the volumetric density PDFs of s = 
ln(p/po) for hydrodynamic runs with different forcing 
(from left to right: solenoidal, mixed, and compressive 
forcing), and different Mach number (from top to bot- 
tom: ~ 3, 5, 10, 20, and 50) at a fixed resolution 
of 256^ grid cells. Each panel shows the time evolution, 
with the initial turbulent state t = 0, with SEE = 0, 
representing the time right before the ffist sink parti- 
cle forms, and when the simulation reached SEE = 5% 
and 20%. The initial turbulent state, t = 0, can be ap- 
proximated with log-normal distributions. Equation 
in the low-Mach number cases, while larger Mach num- 
bers lead to stronger deviations from the perfect log- 
normal shape. This is partly caused by numerical and 
partly caused by physical effects. The numerical effects 
are poor sampling and limited resolution at high Mach 



number ([P rice fc F ederrath|2010| [Price et al.|2011[ jKon- 

standin et al. 2012b|. We explore the resolution depen- 
dence and the influence of different random sampling in 
Appendices [A[ and [C] The physical effect introducing de- 
viations from log-normal statistics is an increased inter- 
mittency for high-Mach number turbulence, producing 
larger skewness and kurtosis in the distributions | Klessen 
2QQQllKritsuket al.|2007||Burkhart et al |2009[ [Federrath 



et al.pOlObl |Hopkins||2Q12bp . Although both numerics 
and physics attect the wmgs of the PDFs, we can draw 
general conclusions from the trends seen in Figure [3j 

At late times, all simulations develop power-law tails 
when the gas collapses to high densities. This effect is 
clearly time-dependent, such that the slope of the tail 
becomes flatter as more and more gas collapses. This 
was seen ffist in smoothe d parti cle hydrodynamics (SPH) 
simulations by Klessen (2000) and la ter confirmed in 



grid si mulations with tracer particles (Federrath et al 



2008a[) and in pure grid simulations ([Pib ^ Burkertl! 
Vazquez-Semadeni et al.'2008MCho fc k im|Ml ||KT] 
et al. 2011a; Ballesteros-Paredes et al. 20111 ICo 



_ _ic 

2012( [Safranek-Shrader et al.| |2012), though with a sig- 
niticantly narrower parameter range than studied here. 
All those prev ious numerical models (except the study by 
Klessen 2000[ which on the other hand had a very low 
resolution of only about 70^ effective SPH particles) were 
limited in that they could not follow the collapse to late 
stages with SFEs much larger than SFE > . be cause 
they did not include sin k particles (see Section ^j j). For 



[Kritsuk 
ins et al. 



instance, 



Kritsuk et al. 



had to stop their simulat 



( ^20 11 a) and Collins et a! 
tions when the first over-c 



2012) 
densities 



started to collapse, because the simulations become pro- 
hibitively slow and the resolution criterion for c ollapsing 
gas is violated at late times (see Section 2.2). Using 
sink particles here, we can follow the development of the 
high-density tail in the PDFs to much later times. To 
guide the eye, we draw power-law PDFs that would be 
obtained for purely radial density distributions. 



p{r) (X r 



(6) 



with = 1, 2, and 3 (dotted, solid, and dashed power- 
law lines in each panel). Here we compare PDF slopes in 
terms of this equivalent radial density profile. This en- 
ables a direct comparison of the slopes in the tails of the 
density PDF from different numerical, theoretical, and 
observational studies, because the slope of the equiva- 
lent radial profile can always be computed from any rep- 
resentation of the density PDF (e.g., volume- weighted or 
mass- weighted, linear or logarithmic density, volume or 
column density). We provide a derivation of the rela- 
tionship between volume and column density PDFs and 
their equivalent radial density profiles below. A compar- 
ison of high-density PDF tails in terms of equivalent ra- 
dial density profiles is furthermore motivated by the fact 
that the dense parts of molecular clouds eventually de- 
velop self-gravitating cores, which can be approximated 
with such radial power-law profiles as studied in theory 
(e.R.,[L arson"1969! 'Penston 1969| |Shu||1977[ iWhitworth 



|& Summers 1985), numeric al simulations (e.g., IFost er 
Chevalier 1993; Ba nerjee fcP udritz 2007; Federrath e t a. 
|2011c|), and obse rvations (jK onyves et al. 2010; Arzouma- 
|nian et al.|[2QTT| iHih et alT 2011; Schneider et al. 201^ 
Here we use volume- weighted PDFs of the natural loga- 
rithm of density, Ps with s = \n{p/ po), as in Equation (^. 
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^-20, sol, ^0=0, 256^ ^=0^ff- 
SFE= 0.0% ■ 

SFE= 5.0% 
SFE = 20.0% 




_ ^-50, sol, 5o=0, 256-' t=Ot^^ 
SFE= 0.0% 
SFE= 5.0% 
SFE=1 9.7% 





^-20 
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Figure 3. Probability distribution functions (PDFs) of the logarithmic density s = \n(p/po) for hydrodynamic (Bq = 0) models with 
solenoidal (left), mixed (middle), and compressive forcing (right) and sonic Mach numbers Ai ~ 3, 5, 10, 20, 50 (from top to bottom) 
at different times: t = 0: state of fully developed turbulence (solid), SFE = 0: right before first sink particle forms (dotted), SFE ~ 5% 
(dashed), and SFE ~ 20% (dot-dashed). The standard deviation of the PDFs increases with increasing Mach number and more compressive 
forcing. All PDFs develop power-law tails at high densities of flattening slope with increasing SFE. The three power-law lines below the 
legend in each panel show slopes of power-law PDFs equivalent to PDFs obtained for volumetric, radial-density profiles p(r) oc with 
K = 1 (dotted), K = 2 (solid), and = 3 (dashed), according to Equation (|9|. 



The PDF Ps is related to the PDF of the density, Pp by 
(e.g., |Li et al. [20031 [Federrath eral]|2011c| ) 



with the volume V oc r^. Substituting this into Equa- 
tion ([7|), we find 



Ps PPp - 



(7) 



Ps oc p ^/'^ (X ex.p{—3s/K,) . 



(9) 



The PDF Pp in turn is a measure of the volume fraction In analogy, we can derive the power-law scaling of the 
with gas of density p, and is thus defined as PDF of column density S oc pr oc r~'*+-^, 



Pp (X dV/ dp = ( dV/ dr)( dr/ dp) cx (g) 



Ps cx dA/dS cx s(i+'=)/(i-'=) 



(10) 
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Figure 4. Same as Figurejs] but for the logarithmic column density 77 = ln(E/Eo) with the mean column density Eq. Error bars indicate 
the uncertainty in the PDFs when three different projection directions (along the x, y, or z axes) of the same cloud are compared. The 
column density PDFs are more noisy, but follow the same trends as the volurnetric density PDFs in Figure [ 3 | However, the standard 
deviation is smaller in prj than ps, because of projection (Fischera & Dop ita)2004[ [Federrath et al.[2010b| [Brunt "et al.|2010b|a[[Seon|2012| ), 
and the power-law tails have a different slope (compare Equations [y] and ^ 



with the area A oc r'^. The logarithmic column density 
PDF of 7^ = ln{T^/T^o), defined in analogy to s, is 

Pr, = (X S^/^i-'^) (X exp [2r]/{l - i^)] , (11) 

relating the slope of the column density PDF to the slope 
K of the equivalent radial density profile given by Equa- 
tion (|6|. Note that the column density PDF slope di- 
verges for = 1, because the column density, S ex r~'^+^ 
is constant in that case, which means that and are 
delta functions. 



4.2.2. Column Density PDFs 

Figure [4] shows the same as Figure |3j but for the log- 
arithmic column density variable rj = ln(I]/So), where 
So is the mean column density of the region considered. 
The power-law slopes for radial, volumetric density pro- 
files are given in analogy to Figure [3] for = 1, 2, and 
3 in each panel. The transformation of slopes between 
the volumetric and the column density PDFs is obtained 
by comparing Equations ([9| and (11). Intermittency and 
sampling effects are more pronounced in the column den- 
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Figure 5. Top row: same as Figure [3] but for models with mixed-forcing , 
1 /iG (middle), and 3 /iG (right) at 512^ resolution. Bottom row: same as to 



at ~ 10 and varying initial magnetic field Bq = Ofleft), 
^ ^ , top row, but for the column density PDFs as in Figure |4] The 

density variance decreases with increasing magnetic field strength, in agreement with the theoretical prediction given by Equation ([5|. 
High-density power-law tails develop in all models, but no clear trend with magnetic field strength is seen. 



sity PDFs than in the volumetric density PDFs. We also 
added error bars to indicate the variation of the r^-PDF 
when three different projections of the same cloud are 
considered. The general trends with different models are 
similar to the volumetric density PDFs of Figure The 
statistics at low column density is more strongly affected 
by sampling and projection effects, but the power-law 
tails at high densities are significant, and can be mea- 



sured in molecular cloud ob servations (e.g., Kainulainen 
et al.||2QQ9| [Schneider etsTl mU). 



stop the simulation as soon as the first over-density went 
into collapse. We find a similar global slope in our model 
with mixed forcing at ^ 5 when SFE ~ 0. However, 
our results for different forcing and sonic Mach numbers 
suggest that the slope of the high-density tail of the PDF 
is not universal. Moreover, the PDF slope is clearly time- 
dependent and flattens in all our models, consistent with 
previous st udies probing a more limited physical param- 



eter range (| Cho fc Kim]|2Qll[ [Collins et al.||2Q12| ). 



TEe power-law slope of the high-density PDF tails de- 
pends on the range chosen to fit and on the SFE. It 
may also depend on physical parameters (different forc- 
ing and/or Mach number, etc.), but — unlike the standard 
deviation of the PDF, Equation ([5| — the power-law slope 
does not seem to depend strongly on the forcing or the 
sonic Mach number of the turbulence. Here, we only 
highlight general trends without fitting slopes, because 
of limited resolution, sampling, and general uncertainties 
in setting the fit range. However, we clearly find that the 
high-density tail in both volumetric and column density 
PDFs of Figures [3] and [4] flattens with increasing SFE, 
or equivalently, the corresponding radial density profile. 
Equation |6), steepens from k, < 1.5 for SFE ~ to 
- 1.5-2.5between SFE - and SFE - 20% in ah nu- 
merical models. Our estimates of agree well with the 
range of equivalent radial power-law slop es inferred in 
recent Herschel observations, ~ 1.5-2.5 (Arzoumanian 
[eral.[[2011[), suggesting SFEs in the range ^FE 

— 



4.2.3. The PDF of Magnetized Clouds 

Figure [s] shows a comparison of volumetric (top) and 
column density PDFs (bottom) for simulations with 
magnetic fields, Bq = (left), laG (middle), and 3 /iG 
(right). According to Equation (pi), the density variance 
should decrease with increasing magnetic pressure. In- 
deed, the PDFs become narrower with decreasing plasma 
or d ecreasing Ma (see Table [T|, as observed in Collins 
et al.[ |2Q12), and quantified in Molina et al. (2012]). The 



high-density power-law slope, however, does not seem to 
exhib i t a cle ar trend with magnetic field strength. Collins 



Kritsuk et al. (2011a) provide a very detailed measure- 
ment of the slope in their single model with mixed forcing 
at ~ 6. They find a global PDF slope corresponding 
to K ^ l.S at the latest time available in their simula- 
ti on. Steeper global pow er-law slopes were not obtained 



et al. (2012) report a weak steepening of the tail with de- 
creasing^ but that depends again on the time and fit 
range chosen. 

Correlating model parameters with measured PDF 
slopes for diflFerent SFE would be a natural way to pro- 
ceed in follow-up studies, in order to learn about the 
state of a cloud by measuring its column density PDF. 
For instance, the standard deviation (width) of the PDF 
shows a clear dependence on the forcing, the Mach num- 
ber, and the magnetic field strength. These dependences 
have been studied theoretic ally and numerically for volu- 
metric den sity PDFs dP assot & Vaz quez- Semadeni|1998| 



Nordlund & Padoan[[ 



EiliKiaiiaBBf^ 



Kritsuk et al. (2011a), probably because they had to 



|et al.||2011 ; Molina et a l.||2012( ^Konstandin et al. |[2012lDf 
and were extended to the study of column densityTTJFs 



The Star Formation Efficiency 



11 



iFisch era^ 

Burkhart 







2012 




that t 




as A4 is increased and as the forcing becomes more com- 
pressive, while we see in Figure [5] that the standard devi- 
ation decreases with increasing magnetic field strength. 
Thus, a promising goal would be to relate PDF parame- 
ters with cloud properties, in order to constrain the evo- 
lutionary stage and physical paramete rs of a cl oud, e.g., 
star- forming versus quiescent (Kainul ainen et al. 2009), 
pressure-confined or not (Kainulainen et al. 2011), in- 
fiuence d by ionizing radiatio n, or other environmental 



effects ( Schneider et al.|2012 ), by measuring column den- 
sity PDl^ s and comparing them with simulations. 

5. VELOCITY, DENSITY, AND COLUMN DENSITY 
SCALINGS OF MOLECULAR CLOUDS 

In this section, we study the spatial scaling of turbu- 
lent velocity, density, and column density. In particular, 
we find that the column density power spectrum is a 
powerful tool to distinguish star-forming from non-star- 
forming clouds. 

5.1. Definition of the Fourier Power Spectrum 

Fourier power spectra are commonly used in turbu- 
lence analysis to study the scaling and correlation of 
a given physical quantity q{£) (e.g., velocity, density, 
or combinations) with spatial scale £. Since the anal- 
ysis is carried out in Fourier space, the spatial scale 
simply transforms to a wavenumber scale k = 27r/i. 
The D-dimensional Fourier transform of q{i) with i = 
is defined as 



q{k) = 



(27rL)^/2 



ik£ aDc 



where we denote the Fourier transform of q{i) with q{k)^ 
where k e 3?^. The spatial scale £ is in the range 
[0, L] with the maximum spatial scale L corresponding 
to the smallest wavenumber 27r/L, i.e., the largest scale 
of a cloud or in a column density map. Note that in this 
definition of the D-dimensional Fourier transform, q has 
units of ^ [qk-^/'^], such that Jq-q'd^k al- 

ways yields a quantity with units of the original physical 
variable q'^. 

With the definition of the Fourier transform. Equa- 
tion (12), the 3D {D = 3) Fourier power spectrum of q 
is given by 



(13) 



as an average of q- q^ over a spherical shell with radius 
k = \k\ and thickness dk in Fourier space, where q^ 
denotes the complex conjugate of q. In analogy, the two- 
dimensional (2D, D = 2) power spectrum. 



P2B{q,k) = {q-q^27Tk)k, 



(14) 



is defined as an average over a ring with radius k and 
thickness dk. As a result of the definition of the 
I^-dimensional Fourier transform, both PsBiQ^k) and 
^2d(^, ^) have the same units, [q'^/k]. 



5.2. Definition of the A-variance Spectrum 

The A-variance technique provides a complementary 
method for measuring the scaling exponent of Fourier 
spectra in the physical domain by applying a wavelet 
transformation ( Stutzki et al.|1998 ). W e here use the A- 
variance tool developed and provided by |Ossenkopf et alT 
(2008) for 2D maps. This tool impl ements an improved 
version of the original A-variance (Stutzki et al. 1998 
Bensch et al. 2001) , capable of treating non-periodic 
edges and takmg mto account variations in signal-to- 
noise, which is important for applying the met hod to 
observational maps (e.g., Schneider et al.|[2011 ). The 
A-variance measures the amount of structure m an ob- 
servational map, by filtering the dataset q{i) with a sym- 
metric up-down-function 0^ (typically a French-hat or 
Mexican-hat filter) of size ^, and computing the variance 
of the filtered dataset as a function of filter size £. The 
A-variance is defined as 



(15) 



where the average is computed over all data points at 
positions x = {xi^X2} for a 2D map (^(x). The operator 
* denotes a convolution, which in practice, is carried out 
in Fourier space. We use the original French- hat filter 
with a diameter ratio of 3.0 as in previous studies using 
the A-variance technique (e.g., Stutzki et aL]|1998| Mac 



'Low Ossenkopf 2000VUssenkopl et al.^ 2001||r)ssenkopl 

& Mac Low 2002; Ossenkopf et al. 2006)"There is some 
uncertainty introduced by the filter type and diameter 
chosen, but for turbulent systems, the A-variance does 
not s ignificantly depend on that choice (Ossenkopf et al. 
2008 ). Here, we apply the A-variance only to the col- 
umn density, i.e., q = in Equation (15). The relation 

2 



between cr^(i;,£) and the Fourier power spectrum of the 



(12) density, P3d(p, ^), is explained in the next section. 



5.3. Relation between Psb cind cf\ 

Assuming a power-law scaling of the 3D density power 
spectrum, 

PsB{p,k)o^k'' (16) 

and a power-law scaling of the 2D column density power 
spectrum 

P2D(S,fc)ocfc^, (17) 



we see from the definitions ( |13[ ) and (14) that they are 
directly related by 



3D 



2kP2B (X /C" (X . 



(18) 



Moreover, considering the spatial scaling of the column 
density S in real space (instead of wavenumber space), 
we have to recall that P2d(S) oc dS^/ dk. Thus, 

(X P2Bk (X Psb oc k^^^ ex r (^+^) ex . (19) 

This result shows that the A-variance power spectrum 
of S defined in Equation ( 15 ) is simply related to the 3D 
Fourier power spectrum oTp by 



Cri(S,^) (xP3d(p,/c) CX CXr" 



(20) 



for PsBip^k) CX k^. Thus, the A-variance slope of S is 
simply the negative slope of the power spectrum of p, if 
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both follow power-law scalings. We will make use of this 
result in Section [71 below. 

5.4. Velocity Spectra 

We start by considering velocity spectra in order to 
compare them with previous results for supersonic tur- 
bulence, and to study their Mach-number and forcing 
dependences. Figure [6] shows power spectra of the ve- 
locity vector V normalized by the sound speed Cg, i.e., 
the local Mach number. Thus, we set g = v/cg in Equa- 
tion (|13|) and define the shortcut v' = v/cs . Since the 



spectra {t = 0, solid lines) a re roughly consiste nt with 



simulations cover a wide range of RMS Mach numbers 
M. ~ 3-50, and hence the absolute power of each spec- 
trum varies depending on Al, we chose to normalize all 
velocity power spectra by their individual integral in Fig- 
ure [6j such that the shape of the spectra (and in partic- 
ular their slopes) can be more easily compared between 
the different numerical models. Moreover, the spectra 
are compensated by /c^, such that a horizontal line would 
correspond to a power-law spectrum ex k~'^ as in Burgers 
(jl948J turbulence. For comparison, we plot two dotted 
lines with p ower-law scalin g ex k~^^^ (Kolmogorov |1941| ) 
and cx k~^ (Burgers"l948). 

Before starting to analyze the dependence of the power 
spectra on the forcing, sonic Mach number and magnetic 
field, we first have to identify the scales on which the 
spectra are robust, numerically converged, and can be 
directly compared. At small wavenumbers, k = 1-3, the 
turbulent forcing directly affects the spectra, bec ause ki- 
netic energy is injected there (see Section 2.1), which 
then self-consistently cascades down to snialler scales 
(larger wavenumbers). At the high- wave numb er end, nu- 



merical resolution and t he bottleneck effect (Falkovich 
1994| Dobler et al.|2003 ) affect the spectra. As shown in 



previous studies, scales smaller than about 30 grid cells, 
i.e., k > fcdiff with k^m ^ 512/30 - 17 for the 512^- 
resolution spectra analyzed h ere, are typical ly affected 



by numerical diffusion (.Kitsionas et al. 
et al.||2010b| [Price Federratli||!i(jloT 



2009 Federrath, 



IHur et al.[[:^Olor 



Federrath et al.| |2011c| [Kritsuk et al.||2011bp . Moreover, 
the bottleneck ettect could in principle attect the spec- 
tra at slightly smaller wavenumbers than /cdiff- However, 
previous studies indicate that the bottleneck effect has 
a minor impact on velocity spectra of supersonic tur- 
bulence (unlike in mildly compressi ble or incompressible 
turbulence; see Dobler et al. 2003). Compressive forc- 
ing, in particular, does not seem to produce a significant 
bottleneck effect at all (see the 1024^ resolution spectra 
of ~ 5 .5 turbulence with sole noidal and compressive 
forcing in Federrath et al.| 2010b). This is likely because 
steeper spectra (e.g., oc k~"^ as in highly compressible, su- 
personic turbulence) have a smaller correction due to the 
bottleneck effect than shallower spectra (e.g., cx k~^^^ as 
in incompressible turbulence) (Falkovich 1994), and so 
the bottleneck effect is weaker m supersonic than in in- 
compressible turbulence. Thus, both forcing {k < 3) and 
primarily numerical diffusion {k > 17) limit the scaling 
range in which all power spectra can be reasonably ana- 
lyzed and compared. The remaining, intermediate range 
is quite small, but still allows us to draw general conclu- 
sions. In order to leave some tolerance on both ends, we 
define the approximate inert ial range here as 5 < /c < 15. 
In Figure [6) we see that the purely turbulent velocity 



power laws, Psb cx k for Kolmogorov| ( 194l| turbu- 
lence below the sonic scale {k > /Cgj, and are generally 
consistent with Psd oc A:~^ on scales larger than the sonic 



scale ( k < k^ ). W e calculate the sonic scale as in [Schmidt 



et al.[ ( |2009[ ) and [Federrath et al.| ( 2010b| ) by evaluating 
the implicit integral definition Psd {v/c^^k)dk = 1. 
This definition states that the kinetic energy {v'^)/2 on 
all scales k > k^ must be equal to the sound energy Cg/2 
(note that here we did not include a factor of 1/2 in the 
definition of the power spectrum. Equation [Tsj so it sim- 
ply cancels on both sides of the integral definition of k^). 
The sonic scale thus separates supersonic scales {k < k^) 
from subsonic scales {k > k^). We indicate k^ by the ver- 
tical lines on top of the abscissas in Figure |6] It shifts to 
smaller scales (higher k) for models with increasing 
which simply means that a larger fraction of resolved 
wavenumbers remains in the supersonic regime. Hence, 
in the ~ 50 runs, the turbulence is supersonic on all 
scales, and the sonic wavenumber is larger than the cut- 
off wavenumber given by the numerical resolution. On 
the other hand, in the A4 ^ 3 runs, the sonic scale is in 
the resolved range of scales, where the spectra are con- 
sistent with cx k~^^^ scaling for k > k^. Even though the 
spectra are quite noisy (because we only consider one re- 
alization of the turbulence at a fixed time) and the range 
of scales for analysis is limited as explained above, we 
find a flattening of the spectra from P cx k~'^ in the su- 
personic regime {k < k^) to P cx k~^^^ in the subsonic 
regime {k > k^). 

Compressive forcing exhibits slightly steeper spectra 
than solenoidal forcing for any given Mach number. This 
is because compressive forcing produces stronger s hocks , 
and hence we expect a scaling closer to Burgers] (1948) 

-2 



scaling, P cx k~^. A slight steepening of velocity spectra 
was previously seen only in simulations with A4 ~ 5-6, in 
which solenoidal forcing yielded a slope of —1.86 =b 0.05 
and co mpressive f orcing a slope of —1.94 ± 0.05 (Fed-[ 
errath et al. 20'09l |20T0b|. The reason for the slightly 



steeper velocity spectra is that compressive forcing ex- 
cites aJ^xge£fe£ti£nof_ooim modes at any given 
M ( Federrath et al.|2011b| ), which — in the limit of a 
hypersonically turbulent velocity field, completely filled 
with shocks — is the essence of Burgers (1948) scaling of 
turbulence. 

Observations of molecular clouds support a velocity 
scaling between k~^^^ (Kolmogorov) and k~'^ (Burgers), 
corresponding to cx and v cx ^^/^ respectively, with 
most clouds ex hibiting a scaling clo ser to Burgers tur bu- 
lence (|Larson|1 981; Sol omon et al.[ri9 87; Falgarone et al.[ 
19921 lOssenkopf Mac Low[[!^00!^[Tneyer Brunt ..2 
^ ^ ^ He 



-Duval et al 



201 ip . Hence, both solenoidal and 
compressive forcing produce velocity scalings consistent 
with observations. 

Our velocity spectra do not strongly change when star 
formation proceeds, as can be seen by comparing the 
t = and SFE = 0, 5%, and 20% curves in Figure [e) 
Only the runs with A4 ^ 3 show a significant dependence 
on SFE, because in those models, global collapse starts 
to increase the velocity power on large scales, /c ~ 1, and 
thus reduces the relative power on smaller scales, be- 
cause of the normalization of the spectra. Thus, velocity 
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sol, ^0=0, 512' £=0^ff- 
SFE= 0%- 
^-5/3 SFE= 5%- 
....■---'f^-2 SFE = 20%- 



Jl-^] 0, sol, Bq=0, 512' 
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Figure 6. Normalized velocity spectra of q = v/cg defined in Equation \13\ and compensated by k'^ for hydrodynamic simulations with 
resolution of 512^ grid cells and sonic RMS Mach number ~ 3 (left), ~ 10 (middle), and ~ 50 (right) for solenoidal forcing (top) 
and compressive forcing (bottom). Different line styles indicate different times or SFEs: t = (the time when turbulence is fully developed 
and self-gravity is switched on; solid), SFE = (right before the first sink particle forms; dotted), and SEE = 5% (dashed) and 20% 
(dot-dashed). Straight dotted lines indicate Kolmogorov (oc k~^^^) or Burgers (oc /c~^) scaling of turbulence in the approximate inertial 
range 5 < /c < 15. Estimates of the sonic wave number ks are shown as vertical lines on top of the abscissa. The sonic scale separates the 
range of supersonic turbulence (k < ks) from subsonic turbulence (k > ks). 



spectra do not seem to be a suitable indicator for star for- 
mation, unles s the collap se is globa l, which is typically 



not observed (fZuckerman fc Palmerl |1974| Zucke rman 



Evans||l974| [Krumholz <§rTan||2QQ7[|Evans et aLn2QQ 9 f\ 



Our results provide an explanation tor the insensitiv- 
ity of the Principal Component Analysis of velocity fluc- 
tuations measured with CO lines in the comparison of 
Rosette (a star- forming cloud ) and G 216-2.5 (a quiescent 
cloud) found in iHeyer et al. (2006). CO traces primar- 
ily the large-scale structure of molecular clouds, because 
the 1-0 and 2-1 rotational transition lines become op- 
tically thick at relatively low column density for both 
^^CO and ^^CO. The results in Figure [g] indeed indicate 
that large-scale velocity analyses of a molecular cloud are 
not particularly sensitive to star formation. In the next 
section, however, we show that the density and column 
density power spectra do depend on star formation, even 
on relatively large scales, and correlate with the SFE. 

5.5. Density and Column Density Spectra 

For the density and column density spectra, we define 
p' = p/ pQ — 1 and = S/Sq — 1 with the mean density 
po and the mean column density Sq. We then set q = p' 
in Equation (13), which defines P3D(p^ ^), and set q = T^' 
in Equation (T3), defining P2D(S^ k). The normalization 
by the mean values is useful, because it enables a direct 



^ See, however, the discussion in |Vazquez-Semadeni et aT] (|2010| 
suggesting that global collapse drives tu rbulence and star torma- 
tion. For in stance, Klessen & Hennebelle (2010) have argued and 
[Federrath et al. (2011c) have shown that collapse and accretion do 
drive turbulence, but at the same time, the turbulence produced 
in this conversion from gravitational energy cannot hold the col- 
lapse. So global collapse would still be observed, unless the stars 
form much more quickly than global contraction occurs, adding 
an additional source of turbulence, possibly capable of holding the 
global collapse or even dispersing the cloud. 



comparison of the spectra, even when the absolute den- 
sity and column density scales are different in the various 
simulation models (see Table [l|. 

Figure [7| shows the same asT'igurejG) but for the den- 
sity and column density fluctuation spectra. Note that 
we expect P3D(p^ ~ 2/cP2D(S^ ^) from Equation (18), 
which is well confirmed in FigurejTl showing 2/cP2d(^ S 
as thick lines and P3D(p^ k) as tnm lines superimposed. 
The match is quite close, as expected for nearly isotropic, 
turbulent density distributions. The gray error bars in- 
dicate variations of P2D(^^^) with different projection 
directions along the x, and z axes, showing that there 
is no significant variation of the column density spectra 
with projection direction. 

We plot density spectra for different times and SFEs, 
t = and SFE = 0, 5%, and 20% in Figure ^ We 
see that the dependence of Ps^ip^k) and P2d(S ,^) on 
SFE is very strong. All spectra have a negative slope a 
in the purely turbulent regime before self-gravity is con- 
sidered in the simulations {t = 0) until the first sink 
particle forms (SFE = 0). We here quote the slope 
a of the 3D density spectrum by fitting a power law 
within the approximate inertial range, 5 < < 15, with 
P3D(p^^) oc k^. Note that the column density power 
spectrum P2d(S',/c) cx k^~^ (see Equation 18), which 
is why we plot 2A:P2d, also exhibiting the slope a. For 
SFE = 5%, this slope becomes positive and increases 
further for SFE = 20% in all numerical models, irre- 
spective of the forcing and sonic Mach number. This 
result and our measurements of the slope a do not de- 
pend on numerical resolution or statistical sampling of 



the turbulence as demonstrated in Figures [13] and [14 
the Appendix. 

The difference of a between t = and SFE = 20% is 
largest in the ~ 3 runs, intermediate for ~ 10, 
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Figure 7. Same as Figure Im but for the column density fluctuation spectra P2D(^^ with q = S/Eq — 1 in Equation \1A\ and multiphed 
by 2k (thick hues) to transiorm the column density spectra to volumetric density spectra P3D(p^ k) with q = p/ po — i m Equation (|13| 
(thin lines). Both 2D and 3D spectra agree as expected from Equation JTsJ, within the uncertainties introduced by different projection 
directions along the y, or z axes, indicated by error bars. Power-law flts m th e approximate inertial range, 5 < /c < 15, show that the 
density spectra evolve from negative power-law slope a deflned in Equation ( |16| to a positive slope, when star formation proceeds in the 
clouds (SEE > 0), irrespective of the particular model parameters. This makes the density and column density spectra a powerful tool to 
distinguish the evolutionary state of a cloud and to estimate the SEE in observations, as will be discussed in Sections [6] and [T] below. 
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Figure 8. Same as Figure[7| but for ~ 10 models with mixed forcing and varying magnetic fleld strength Bq = ^ (left), 1 pG (middle), 
and 3 pG (right). The slope a of the density spectra changes from negative to positive for star-forming clouds (SEE > 0), irrespective of 
the magnetic fleld. 



and smallest in the ~ 50 simulations. Comparing dif- 
ferent forcing, we see that on average, solenoidal forcing 
produces systematically larger slopes a than compressive 
forcing at any ^ive n Mach number. This has been no- 
ticed previously by Federrath et al. (2009), but only for 
M. ~ 5-6 simulations in the purely turbulent regime, and 
is related to the different fractal dimensions obtained 
for solenoidal forcing {D^ ~ 2.6) and compressive forcing 
{Df ~ 2.3). Here we see that the higher the Mach num- 
ber, the smaller are the differences between solenoidal 
and compressive forcing, which is plausible, because we 
expect both forcings to produce a statistically identical 
network of pure shocks for oo, which is the the- 

oretical limit of Burgers turbulence. We also see that 
in this limit, the differences in slope with time and SFE 
become smaller, because gas at very high Mach numbers 
is so strongly compressed in shocks and dense filaments, 
that collapse does not change the density structure as 
much as in the case of lower Mach numbers. However, 
the fact that the slope a is negative for SFE ~ and 
becomes positive for SFE ^ is robust, even when the 



Mach number is increased from ~ 3 to 50, and even 
when the forcing is varied from one extreme (solenoidal 
forcing) to the opposite extreme (compressive forcing). 

In Figure |8j we see the same general trend of increasing 
slope a with SFE in the MHD simulations with A4 ~ 
10 and mixed forcing. The density power spectra do 
not seem to depend systematically on the magne tic field 
strength, consistent with the conclusion drawn in |Collins| 
et al. (2012). The spectra for B = 1 jaG and 3 jaG have a 



slope a that is slightly larger (by about 0.2 on average) 
than in the B = case, which is a relatively small, but 
noticeable difference. However, the change in slope from 
a < for SFE - to a > for SFE > is clearly seen 
in the MHD models in Figure |8j as before in Figure [7| for 
different forcing and sonic Mach number. 

All density spectra in Figures [7| and |8] are broadly con- 
sistent with previous studies exploring a limite d subset 
of the parameter space analyzed in this study (Beresnyak 



oi tne parameter space analyzed m tnis s tuay ( i^e resnyaK 
etaI72005! 'Kowal et al.'2007^ "Kritsuk et al. 2007} |l^'eder-' 
rath et al. 2009; Collins et al. 2012). Here we show that 
all density spectra exhibit a Battening with increasing 
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Mach number for both solenoidal and compressive forc- 
ing, con sistent with a previ ous study of purely solenoidal 
forcing ( |Kim fc Ryu||2QQ5 ). However, the most striking 
result is that all density spectra exhibit negative power- 
law slopes for SFE ~ 0, but turn to positive slopes for 
SFE > 0, making them a potentially powerful tool to 
distinguish star-forming clouds (SFE > 0) from clouds 
that did not (yet) form stars (SFE = 0). Moreover, the 
slope of the density spectra increases roughly monoton- 
ically with increasing SFE (explored in more detail be- 
low). Here we show that the switching in slope from 
negative to positive a of the density and column density 
spectra is universal, holding in clouds with A4 = 3-50 for 
solenoidal and compressive forcing, as well as for different 
magnetic field strengths. This universal behavior of the 
density spectrum is a signpost of gravitational collapse 
when star formation proceeds. As more and more power 
in the density field moves to small scales where star for- 
mation occurs, the density spectrum rises on those scales 
and causes the significant change of the slope a that we 
see here. 

5.6. Delta- variance Spectra of the Column Density 

As a more practical approach to measure the density 
power spectrum in molecular cloud observation s, we now 
study the A-variance spectra, defined in Section [5^ The 
A-variance tool provided by Ossenkopf et alj ( |2QQ8 ) can 
treat spatial variations in the signal-to-noise ratio, as well 
as non-periodic boundaries, making it a useful tool to 
measure the column density spectrum in observational 
maps. For pure power-law distributions, the A-variance 
of the column density, a 'iiTi.i) fflves informa tion equiv- 
alent to Fourier spectra ( Stutzki et al.||1998 ), and is re- 
lated to the 3D density power spectrum Ps^ip^k) via 
Equation ^Q] ) . 

Figure [9] shows the same models and times as in Fig- 
ure [7| but for the A-variance spectra of the column den- 
sity contrast (T^(E/I]o, ^), where we divide by the mean 
column density to facilitate the comparison between 
different models. As for the Fourier spectra of the den- 
sity and column density in Figurej7| we see that the slope 
of <j^(S/So,^) depends sensitively on local collapse and 
star formation. A similar trend was repor ted in simu- 
lations analyzed by Ossenkopf et al.| (2001), yet with a 
much more limited dynamic range and set of parame- 
ters. The basic e ffect was already seen in that study, but 
[Ossenkopf et al. ([2001 ) did not relate the slope a to the 
^'E, because ot the relatively large uncertainties in their 
simulations. For the present set of simulations, however, 
we can determine the slopes a reliably in the approxi- 
mate inert ial range, 0.06 < £/L < 0.2, without being 
affected by numerical resolution issues (see Figure 13 in 
the Appendix). 

We recast the slope of the A-variance into the slope 
a of PsD (p, k) via Equation ([20| to make them directly 
comparable to the Fourier spectrum slopes. Comparing 
the slopes a in Figures [7| and [9) we find good agreement. 
The difference in slopes between the A-variance and the 
Fourier spectra is typically less than ±0.3, caused by 
statistical uncertainties in the fit range and by the fact 
that the inertial range is only approximately described 
by a power law. However, the changes in slope with 
increasing SFE are clearly significant and consistent with 
the results obtained in Figure [7| 



Moreover, we find that the A-variance slopes of runs 
with different magnetic field also agree well with the 
Fourier spectra shown in Figure [S] The general conclu- 
sion from all spectra is thus robust: the slope a switches 
sign from negative to positive when star formation pro- 
ceeds (SFE > 0), irrespective of the sonic Mach number, 
the forcing, or the magnetic field strength. 

6. ESTIMATING THE SFE 

In the previous section, we showed that the spectral 
slope a of all density power spectra switches from nega- 
tive to positive when star formation proceeds. We now 
study the quantitative dependence of a on SFE. Fig- 
ure flT] shows (a) averaged between the Fourier spectra 
(Figures [t] and |8| and the A-variance spectra (Figures [9] 
and P!q|) as a function of SFE. We see that a increases 
with^E from a < for SFE - to a > for SFE > 0, 
except for the ~ 50 models (top panel), where a ^ 
for SFE ~ 0. These high-Mach number models also show 
the slowest increase, from a ~ for SFE ~ to a ~ 0.8 
for SFE = 20%. The M^3 models, on the other hand, 
show a much stronger difference between SFE ~ and 
SFE > 0. This can be explained by the fact that the 
density structure changes significantly during collapse in 
the low-Mach number models, while the global structure 
is almost unaffected during star formation in the high- 
Mach number cases, because the latter already contain 
most of the mass in high-density filaments caused by tur- 
bulent compression (see Figures [l]and[2|. 

For a given sonic Mach number, we see that solenoidal 
forcing always produces higher a than compressive forc- 
ing (top panel in Figure 11). However, the differences 
between the forcings become smaller as the Mach num- 
ber increases. This is quite plausible, because in the 
hypersonic limit {M oo), we expect both solenoidal 
and compressive forcing to produce maximum compres- 



sion, r esulting in a network of pure shocks as in Burgers 



( |1948[ ) turbulence. The differences between solenoidal 
and compressive forcing with ~ 50 are indeed not sig- 
nificant, indicated by the overlapping error bars in those 
two cases (circles and stars). Real molecular clouds have 
sonic Mach numbers in the range A4 ~ 2-50, with most 
clouds exhibiting ~ 10, where differences between the 
forcings become pronounced when SFE > 5% (triangles 
and squares). For Mach numbers A4 ^ 3 (crosses and 
diamonds), the differences between solenoidal and com- 
pressive forcing are significant in both the purely turbu- 
lent (SFE ~ 0) and star-forming regimes (SFE > 0). 

The bottom panel of Figure 11 shows that adding a 
magnetic field does not change a significantly, except 
in the late stages of star formation when SFE > 10%. 
There, the higher the magnetic field, the slightly larger 
becomes a. 

In order to better distinguish models with different 
physical parameters and to provide a quantitative mea- 
sure of the relation between a and SFE, we made empir- 
ical fits to the data in Figure [TT] We use the following 
model function for fitting. 



(21) 



a = po - exp {pi - P2 X SFE[% 



with parameters Po, Pi, and p2, because of its simplicity 
and because it can be inverted. The positive-definite 
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Figure 9. Same as Figure[7| but for the A-variance spectra of the column density contrast, cr^(E/Eo,^), defined in Equation ( |l5| >. Note 
that the slopes a quoted in the legend were measured in the approximate inertial range, 0.06 < ^/I/ < 0.2, and recast into the slope of the 
3D density Fourier spectrum, P3d(P: k) oc k^. The slopes computed from the A-variance spectra are similar to the slopes in Figure [t] as 
expected for approximate power laws (see Equation |20|). 



10" 
100 
1 

1 

0.1 

0.01 

0.001 
10"^ 



Jlr^] 0, mix, 5n=0, 51 2-' 



SFE=20%- 
SFE= 5%- 
SFE= 0%.. 
t=Ot,- 



Jl'^'i 0, mix, 



:1/xG, 512^ 




slope 
slope 
slope 

slope 



= + 0.9---" 
= + 0.4- - 

= -0.6 

= -0.6 — 



slope a = + 1 .2 -■- 

slope a = + 0.3 - - 

slope a = -0.6---- 

slope a = -0.6 — 



Jl^] 0, mix, 5o=3/xG, 512^ 




slope a = + 1 .5 -■- 

slope cx = + 0.5 - - 

slope a = -0.4 " 

slope a = -0.6 — 



0.01 



0.1 
t/ L 



1.0 0.01 



0.1 
t/ L 



1.0 0.01 



0.1 
t/ L 



1 .0 



Figure 10. Same as Figurelo] but for 7W ~ 10 models with mixed forcing and varying magnetic field strength Bq = (left), 1 /xG (middle), 
and 3 fiG (right). As for the nydrodynamic models of Figure [9] the MHD models indicate a strong dependence of the slope a on SFE. 



Table 2 

Fit Parameters for Each Model in Figure [TT] 



HD models (Figure 11 top) 



solenoidal, 
compressive, 

solenoidal, 
compressive, 

solenoidal, 
compressive. 



M ' 
M ' 
M ' 
M ' 
M ' 
M ' 



Po 



Pi 



P2 



Bo = /iG 
Bo=Oi2G 

: 0/iG 
: 0/iG 

50, Bo=Oi2G 
50, Bp = 0/iG 



3, 
3, 

10, Bo 
10, Bo 



MHD models (Figure 11 bottom) 



mixed, A4 ' 
mixed, A4 > 
mixed, Ai > 



1.990 
1.288 
1.633 
1.149 
0.988 
3.437 



0.788 
0.653 
0.721 
0.332 
0.014 
1.216 



0.097 
0.117 
0.136 
0.097 
0.099 
0.010 



somewhat higher than ±0.2 when appUed to real obser- 
vations, so any further uncertainties (of systematic or 
statistical nature) must be added accordingly. Sources 
of such uncertainties and limitations of the present ap- 
proach are discussed in Section [Sj The fit parameters po^ 
Pi, and p2 are listed in Table [2] for each numerical model 
in Figure [TT] The fit curves are shown in Figure jT as 
solid lines^jray shaded regions centered on the best-fit 
curves indicate an uncertainty range of ±0.2 in a. 



PO 



PI 



P2 



10, Bo = 0/iG 
10, Bo = li2G 
10, Bo=Si2G 



0.921 
1.175 
1.478 



0.389 0.356 
0.450 0.192 
0.642 0.169 



The dependence of a on SFE seen in Figure [TT] and 
provided by Equation (21 ) with the fit parameters in Ta- 



Notes. The model function for fitting is given by Equation ( |21| . 

inverse function 

SFE[%] = max {0, [pi -ln{po-a± 0.2)] /ps} (22) 

thus provides an estimate of the SFE when a is mea- 
sured from an observational map. We included an un- 
certainty in a ± 0.2, reflecting the typical uncertainty of 
a in the simulations to establish the relation in Equa- 
tion (22). In practice, the uncertainties in a are likely 



ble [2] suggests that by measuring the spectral slope a in 
an observational map of dust or molecular column den- 
sity (e. g., by using the A-variance technique described in 
Section 5.2), one can distinguish star-forming (SFE > 0) 
from quiescent (purely turbulent) clouds (SFE ~ 0). 
Moreover, if an independent estimate of the sonic Mach 
number is available (e.g., Burkhart et al.,„2009J and if 
the forcing can be constrained ([Brunt et al. ^2' 



2010b|c 



Brunt|2 010!lK Qnstandin et al.|2U12b||Kainulainen ^ la] 
2012[ Seon 201^, measuring a and comparing it with tb 
present set of numerical simulations may give a hint on 
the SFE of an observed cloud, by inverting the curves 
shown in Figure 11 This inversion is given by Equa- 
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Table 3 

Estimates of the Slope ol of the Density Power Spectrum 
^3d(P5^) cx: k*^ and Inferred SFEs for Galactic Clouds. 



Figure 11. Spectral slope {o;}, averaged between the Fourier spec- 
tra (Figures[7|andj8} and the A-variance analysis (Figures[9]and[l0} 
as a function of SFE. The to'p "panel shows hydro dynamical models 
with solenoidal and compressive forcing at sonic ~ 3, 10, and 
50, while the bottom panel shows MHD models of mixed forcing at 
A4 ~ 10, but with different magnetic field strength Bq = 0, 1/iG, 
and 3 /xG. Error bars include the uncertainties from three different 
projections (along x, y, and z) and the uncertainties from averag- 
ing between Fourier and A-variance spectra. The vertical dashed 
line at SFE = separates the turbulent regime from the regime 
of star formation (SFE > 0). Overall, the slope a monotonic ally 
increases with SFE in each model. Model fits using Equation (|21| 
with fit parameters listed in Table [2] are shown as solid lines. Gray 
shaded regions indicate deviations irom the fits by o; ± 0.2. 



tion (22) with the fit parameters in Table ^ Similarly, 
if SFE and A4 were known for a given cloua, Figure 11 
could be used to constrain the forcing in that cloud. 

7. COMPARISON WITH OBSERVATIONS 

Figure [TT] showed that the spectral slope a is sensi- 
tive to the^FE. Thus, by measuring that slope in in- 
terstellar cloud regions, we can constrain the SFE of ob- 
served clouds. Table [3] lists estimates of a measured with 
different tracers and for different cloud regions, taken 
from various sources in the literature. We find slopes in 
the range a ~ —1.6 to +1.6, and sort each observation 
m Table [s] by a in ascending order (from top to bot- 
tom). For each of the following observational estimates, 
we include both the uncertainties given by the inversion. 
Equation ( [22| ), combined with the fit uncertainties from 
the observation quoted in Table [s] (see also the discussion 
in Section [s] on other uncertainties and limitations). Our 
observational estimates of the SFE are thus conservative 
and strictly yield only upper or lower limits for the SFE. 

The first entry in Table [3] is from a Fourier analysis of 
the 21 cm line of atomic hydrogen of the high Ga lactic 
latitude Ursa Major 'cirrus' by ,Miville-Deschenes et al. 



Cloud/Region 


Ref. 


Tracer 


a 


SFE 


01) 


Ursa Major cirrus 


(1) 


HI 21cm 


-1.6 ±0.2 


~0% 


02) 


Polaris Flare (IRAM) 


(2) 


12CO 1^0 


-1.3 ±0.2 


-0% 


03) 


Polaris Flare (KOSMA) 


(2) 


12CO 2^1 


-1.1 ±0.2 


-0% 


04) 


Polaris Flare (FCRAO) 


(2) 


i^CO 1^0 


-1.0 ±0.2 


-0% 


05) 


Polaris Flare (IRAM) 


(2) 


12CO 2^1 


-0.9 ±0.2 


<1% 


06) 


Vela 


(3) 


dust extinc. 


-0.9 ±0.2 


<1% 


07) 


Polaris Flare 


(4) 


i^CO 1^0 


-0.8 ±0.1 


<1% 


08) 


Perseus/NGC 1333 


(2) 


13CO 1^0 


-0.8 ±0.2 


<1% 


09) 


Coalsack 


(3) 


dust extinc. 


-0.8 ±0.3 


<2% 


10) 


Orion A 


(2) 


i^CO 1^0 


-0.7 ±0.2 


<2% 


11) 


Orion B 


(2) 


13CO 1^0 


-0.6 ±0.2 


<2% 


12) 


Polaris Flare (CfA) 


(2) 


12CO 1^0 


-0.6 ±0.2 


<2% 


13) 


Lupus 


(3) 


dust extinc. 


-0.6 ±0.2 


<2% 


14) 


Mon OBI 


(3) 


dust extinc. 


-0.6 ±0.2 


<2% 


15) 


Mon OBl/NGC 2264 


(3) 


dust extinc. 


-0.5 ±0.1 


<2% 


16) 


Rosette 


(3) 


dust extinc. 


-0.5 ±0.3 


<4% 


17) 


IC 5146 


(3) 


dust extinc. 


-0.4 ±0.2 


<4% 


18) 


Mon R2 


(2) 


13CO 1^0 


-0.4 ±0.3 


<4% 


19) 


Mon OBl/NGC 2264 


(2) 


i^CO 1^0 


-0.4 ±0.3 


<4% 


20) 


Orion A+B 


(3) 


dust extinc. 


-0.3 ±0.2 


<4% 


21) 


Pipe 


(3) 


dust extinc. 


-0.3 ±0.2 


<4% 


22) 


Corona Australis 


(3) 


dust extinc. 


-0.3 ±0.3 


<5% 


23) 


Perseus 


(3) 


dust extinc. 


-0.2 ±0.2 


<5% 


24) 


Mon R2 


(3) 


dust extinc. 


-0.1 ±0.2 


<6% 


25) 


Taurus 


(3) 


dust extinc. 


0.0 ±0.2 


<10% 


26) 


W3 


(3) 


dust extinc. 


0.0 ±0.2 


<10% 


27) 


Chamaeleon 


(3) 


dust extinc. 


0.0 ±0.2 


<10% 


28) 


Cygnus X 


(3) 


13CO 1^0 


0.2 ±0.2 


<17% 


29) 


Serpens core 


(5) 


dust emiss. 


1.6 ±0.4 


>8% 



Notes. The upper and lower limits of the SFE quoted in the 
last column were estimated based on comparing the range of spec- 
tral slopes a (including the maximum uncertainties in Equation [22] 
combined with the additional fit uncertainties for a from the ob- 
servations; see discussion in Section [sj with the numerical simu- 
lations in Figure |11| covering a wide range of physical conditions 
(A4 ~ 3-50 with solenoidal and compr essi ve forcing, and different 
magnetic field strengths). Equation \22\ together with Table j2] 
provid e di rect relations SFE(q;) fitted to the simulation data m 
Figure [TT] References: (1) Fourier analysis by Miville-Deschenes 
^et_al . (liuUcJ); (2) Reanalysis of A-variance estimates in Bensch et ah 
^(3) Reanalysis of A-variance e stimates in Schneider et al 



, |; (4) A- variance by |Stutzki et al., ( ,1998 ); (5j A-variance by 

lOssenkopf et al.| ( |200l| . 



([20031). They find a spectral slope a - -1.6 ±0.2. Com- 
paring this value, including the uncertainties in both sim- 
ulations (Figure 11 and Equation 22) and observations 
(Table [3|, we estimate that the Ursa Major cirrus has 
an SFE~^ 0. This is because even the maximum slope, 
a ~ —1.4 is significantly negative and any of our numer- 
ical simulations did not produce stars for a < —1. Al- 
though molecular condensations have been detected in 
the Ursa Major cirrus, it appears to have essentially no 
star formation activity, likely beca use of the high viria l 
parameters estimated for this cloud (I de Vries et al.llQST ). 
This is in good agreement with our findings ol 'Sb b] ^ 0, 
based on the comparison of the observed slope a with 
our numerical simulations. 

Entries 2-5, 7, and 12 in Table [3| are for observa- 
tions of the Polaris Flare by Stutzki et al. (1998) and 
Bensch et al. (2001") with different mstruments (IRAM 



KOSMA, J^URAU, CfA) and different molecular tran- 
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sitions of ^^CO and ^^CO. We see some variation of 
the derived slope a ~ —1.3 to —0.6, depending on the 
telescope and tracer used, as well as on the spatial scales 
over which the slope was determined ( [Bens ch et al. 2001 ). 
Thus, we reanalyzed the A- variance data in Bensch et al. 
f2001), fitting average po wer laws over a wide r range of 
scales than considered in Bensch et al. (2001), in order 
to obtain a rough estimate of the slope a, less affected 
by the particular choice of scales for the fit. We find 
that our reanalysis mostly agrees with the estimates of 
a in Bensch et al. (2001). We see that the overall slope 
a in the Polaxis^ Flare is of order —1, from which we es- 
timate SFE ~ based on the model fits of Figure [TT] 
and Equation (22) in the Polaris Flare, with a forrnal 
upper limit of SFE < 2%. This is in good agreement 
with observations of the Polaris Flare, showing no sign 
of star formation activity, although the Saxophone region 
does seem to indicate some trace of gravitational contrac- 
tion (N. Schneider 2012, private communication), which, 
however, has not led to any star formation in the cloud 
yet. 

Entries 6, 8, and 9 are for dust extinction and CO ob- 
servation s of Vela, NGC 1333 in P erseus, and th e Coal- 
sack by ISchneider et al.| ( |201l| and lBensch et al.| (|200l| . 



compute SFE ~ SFR x far /-^r with an assumed star 



based on average fits to the A- variance plo ts in Schnei- 



der et al. 



account a 



2011 



broai 



and Bensch et al. (2001), taking into 
range of scales and all the error bars 
of those analyses. The purpose of our reanalysis is to 
obtain general trends of the slope a, without co nstrain- 
ing the fits to re lativ ely small scales as d one in |SchneT 



der et al.| (12011") and [Bensch et aT] (|200lD. This is why 



we find slightly different slopes here compared to their 
small-scale fits. These values of a indicate SFE ~ with 
an upper limit of SFE < 1% in Vela and NGC 1333, 
and SFE < 2% in Coalsack. The latter is in agreement 
with an independent estimate of star forma tion activity 



in Coalsack , suggesting relatively low SFE ( [Kainulainen 
et al.|[2009D . 



Estimates of a for the following clouds and cloud re- 
gions (entries 10, 11, and 13-24) range from a ~ — 0.7± 
0.2 to -0.1 ± 0.2, suggesting SFEs in the range SFE - 
to a maximum of 6%, when comparing to the range of 
SFEs obtained in the simulations for such a (see Fig- 



ure 



11 , Equation 22, and Table [2|. Those clouds include 
Lupus, Rosette, Pipe, Corona^ustralis, and Perseus, 
which indeed seem to be consistent with low to mod- 
erate star formation ac tivity (compare, e.g.. Table 1 in 
Kainulainen et ar]|2009 ). A notable exception may seem 
Orion A-hB, whicn contain a large number of YSOs, in- 
dicating high star formati on act i vity. I n a recent survey 
of Orion A+B, M egeath et al.| (I2012D indeed detected 
3479 YSOs. However, Orio n A+B are als o quite massive 



with Mc - 1.9 X 10^ M0 ( [Wilson et al.||2005| ). Assum- 
ing that YSOs produce stars with a mean mass of about 
(m) ^0.3 Mq (which is in between the mean masses of 
the present-day stellar mass function, (m) = 0.23 M©, 
and the initi al mass f unction, (m) = 0.37 M©, for disk 
galaxies; see ^Chabrier^200 3 ) , we can estimate the SFE of 
Orion A+B with SFE - 3500 x O.3M0/I.9 x 10^ - 
0.6%. A similar estimate based on YSO counts and as- 
suming a median inass for each YSO of 0.5 is pro- 
vided in Lada et al. (2010). From their Table 2, we can 



formation timescale tsF = 2Myr ([Evans et al.|2009 ; Hei- 
derman et al.[|20l'Q Lada et al. [2010p . Using the values 
quoted lor the low extinction threshold, > 0.1 mag, 
covering almost the whole cloud, we obtain SFE ~ 2.1% 
for Orion A and SFE ~ 0.4% for Orion B, consistent with 
our own estimate above. Those SFEs obtained from YSO 
counts, however, are likely lower limits, because of the 
limited detectability of YSOs on small spatial scales in 
crowed, dense regions. The YSO completeness achieved 
with Spitzer may only be 50% averaged over the inner 
0.5 pc of the Trapezium r efflon in Orion, and only 30% 



within a ra dius of 0.2 pc ( Evans et al.|[2009 Gutermuth 
et al.|[20l7| ). Thus, SFE estimates based on Y ISO counts 



The values of a ^ — 0.9 i 0.2 to —0.8 ±0.3 were estimated (Pineda et al. 



can be quite uncertain. However, they do fall in the range 
SFE = 0%-4% that we obtain for Orion A+B by using 
the slope a listed in Table [3] Thus, even if no informa- 
tion on YSOs is available for a given cloud, we can still 
get a lower and upper limit of the SFE from the slope 
a, which is directly measurable from the column density 
map of a cloud alone. 

Taurus, W3, and Chamaeleon (entries 25-27) all have 
a = 0.0 ± 0.2, indicating SFE = 0%-10%. For instance, 
Taurus has a total molecular mass of about 1.5 x 10^ Mq 



2010 ), and more than 300 YSOs were de - 



tected (see I'able 1 in Kainulainen et al. 2009 Elias 1978 ). 
Similar to our estimate lor Orion A+B above, we hnd 
SFE - 300 X 0.3 M0/I.5 X 10"^ Mc^ ^ 0.6% Another es- 
timate based on YSO counts in Lada et aT] (|2010) yields 
SFE ~ 1.1% (again using their low extinction threshold 
value to define Mc). Both are consistent with our for- 
mal upper limit of SFE < 10% from a, but significantly 
smaller than the mean, SFE ~ 5%. Taurus is indeed 
somewhat special, in that it may be more dominated 
by magnetic fields than other clouds. Recent observa- 
tions suggest that turbulence in Taurus is trans- Alfvenic 
( jHeyer fc Brunt,, 2012 J , leading to lower SFRs than in 
super- Allvenic clouds (see Paper I). The rather low SFE 
in Taurus may thus be a result of the relatively strong 
magnetic fields there. Our estimate of SFE = 0%-10% 
from the slope a is slightly biased toward larger values 
of SFE, because stronger roagnetic fields also produce 
slightly larger a (see Figure [TT] bottom panel). 

Cloud refflon 28 in Tab le |^is a A-variance measure- 
ment by Schneider et al.j ([2011) of the Cygnus X giant 
molecular cloud from a high-resolution ^^CO integrated 
intensity map. Cygnus X is one of the most massive and 
most active star-forming regions in the Gala xy, includ- 
ing both o ld and young populations of stars fSchne iderj 



et al. 2006 ). On well-resolved scales, between 4 and 20 pc, 
we tind a power-law slop e a ^ 0.2 ± 0.2 from Figure 2 
in (Schneider et al. (20TT|. Comparing this slope to our 
sirnulation models m l^'igure[TT| we estimate that Cygnus 
X has an SFE in the range ^E = 0%-17%. The lower 
limit of the SFE is formally obtained from the lower limit 
of a, while we actually know that Cygnus X is forming 
stars, so SFE > 0. In particular, the ma ssive DR21 fila- 
ment in Cygnus X shows signs of infall ( Schneider et al.[ 
2010), which is likely the reason for Cygnus X exhibit- 
ing a slightly positive a on length scales representative 
of that filament (length ^ 20 pc, width ~ 5pc) a nd sub- 
filaments collapsing toward the main filam ent ([Schnei- 
der et al. 2010[ jHennemann et aL 2012). The range 
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SFE = 0%-17% found for Cygnus X is relatively large, 
but — to the best of our knowledge — the only currently 
available estimate of the SFE for the whole Cygnus X 
giant mol ecular cloud coniplex. 

Finally, Ossenkopf et al. (2001) studied t he A- variance 
spectrum o f the Serpens core observ ed by Testi fc Sar- 
gent (1998). Ossenkopf et al. (2001) estimated a slope 
a = 1.6 ±0.4, which is clearly positive and larger than 
any of the previous cloud regions studied, indicating 
SFE > 8% (last entry in Table [3|. Based on the total 
available gas mass and the mass spectrum of condensa- 
tions reported in Testi & Sarg entl (|1998[ ), we estimate an 
independent lower limit o f the~5t E in the Ser pens core 
of SFE > 5%. Recently, [Maury et al.| ( |2011| ) reported 
SFE - 7%-15% in Ser pens South and wkO ^nother in- 
dependent estimate by |01mi & Testi| (|2002|) suggests an 
SFE ~ 25%-50% in the region^here most protostellar 
cores and protostars are located. All independent esti- 
mates agree with the lower limit, SFE > 8%, that we 
find here based on a. It is indeed plausible that the 
Serpens core region has a relatively high SFE, because — 
unlike the large-scale cloud regions studied above — the 
Serpens core is a small-scale, high-density region, where 
eventually we expect a relatively large fraction of the 
gas to be accreted by stars. The Serpens core region is 
close to the limit where the SFE approaches the local 
efficiency e, i.e., the mass fraction of gas of an individ- 
ual dense core that will actually fall onto the protostar. 
Theoretical, numerical, and observational estimates sug- 
gest that e = 0.25-0.7 (Matzner fc McKee 2000j Wang 
et al. 2010; Beuther et al.||2002[). l^or instance, | Wilkin; 



statistical and systematic uncertainties are not included 
in Equation ( 22 ) and must be added accordingly. 



jrCada,(,1983pTina" 
'find 07 



e - ITB^TS and 'Ol mi fc Testij 

25-0.5. Comparing the SFR and gas column 



densities in our numerical simulations with Galactic ob- 



servations by Heiderman et al. ( |2010 ) in Paper I, we 
found best-fit values e = 0.3-0.7, which seem to suggest 
typical local efficiencies of e ^ 0.5, in good agreement 
with the aforementioned independent estimates. Thus, 
about half of the gas in a dense core is accreted by the 
protostar, while the other half remains in the envelope 
or is driven out of the core by winds, jets, and outflows. 

8. UNCERTAINTIES AND LIMITATIONS 

8.1. Relation between a and SFE 

In Section [6] we calibrated a simple, empirical relation 
given by Equation ( p2| ) between the spectral slope a of 
the column density and the SFE in the numerical sim- 
ulations of Section [2j This relation can be applied to 
observational maps in which a can be measured with 
Fourier techniques such as the A-variance method, to 
obtain an estimate of the SFE in the observed region. 



Equation (22) only includes an uncertainty in a ± 0.2 as 
suggested by our calibration with numerical simulations. 
The practical uncertainties in a when measured in an ob- 
servational map, however, can be higher than that. Such 
uncertainties include finite telescope resolution, signal- 
to- noise, and sensitivity limits for a given observational 
tracer of column density. For instance, dust emission, 
extinction, or molecular line emission all suffer from sen- 
sitivity limits in both the high and low column-density 
regimes (e.g., molecular lines become optically thick at 
high column density or detectors are simply not sensi- 
tive enough to trace very low column density). Such 



8.2. Global versus Local Efficiencies 

A primary goal of this paper is to understand the tran- 
sition from the global SFE on large scales (e.g., a whole 
molecular cloud) to the local core formation efficiency 
e. The latter denotes the fraction of gas in a dense core 
that eventually goes into the star. This is not 100%, be- 
cause some fraction, 1 — e, of the gas in the core does 
not fall onto the protostar, but is driven out of the core 
by outflows, winds, and/or jets during the formation of 
the protostar in the center of the accretion disk. Without 
specifying what mechanisms (e.g., radiative and/or mag- 
netic) cause the gas flow out of the core or how it actually 
happens, we simply note that we expec t a local efficiency 
e ^ 0.3-0.7 dMatzner fc McKee|[2QQQl Paper I, and ref- 
erences therein). However, even the best high-resolution 
observations of individual cores are just beginning to in- 
vestigate this local regime of star formation and constrain 
e. In surveys of structure and star formation activity of 
molecular clouds on intermediate and large scales, such 
cores appear as sites of high column density, but their 
internal structure is often not resolved. Thus, the actual 
amount of gas contributing to star formation is uncer- 
tain and so is the actual total SFE. We note that we 
did not take into account the local efficiency e in Equa- 
tion (22), because the simulations did not include any 
form oi feedback such as winds, jets, or outflows from 
the protostellar objects. Such feedback effects are not 
accounted for and hence add to the uncertainty in the 
true SFE when Equation ( 22 ) is applied to observational 
data. 

Given the local efficiencies e ^ 0.3-0.7, it would seem 



reasonable to adjust the inferred SFE via Equation (22) 
down by the factor e to estimate the actual total SEE of 
a cloud. In the application of Equation (22) to the obser- 
vations presented in Section [7| however, we did not make 
such an adjustment, because we think that the same ad- 
justment would be necessary in the case of YSO counts 
to infer the SFE. This is because the effect of the lo- 
cal efficiency e is typically also not accounted for when 
estimating the typical mass of a YSO, since those ob- 
servations are limited in resolution and likely include a 
significant fraction of gas/dust in the core that may not 
all be accreted by the star. Thus, the local efficiency 
e due to feedback introduces uncertainties as discussed 
above, which, however, similarly apply to SFE estimates 
obtained with YSO counts. 

8.3. Definition of the SFE in Closed and Open Systems 

We furthermore note that the simulations studied here 
are highly idealized and can suffer from resolution is- 
sues. However, the effects of limited resolution are rela- 
tively small for our measurement of a and we do under- 
stand them quite well (see Appendix |B|). A more severe 
issue is the choice of boundary conditions in the simu- 
lations. They are periodic in all directions. We discuss 
the limitations of this and related, different approaches 
to treat boundaries in Paper I. For the present study of 
the SFE, the main limitation lies in the fact that our 
computational boxes are closed, so there is no gas in- 
flow or outflow. Real molecular clouds are not isolated. 
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closed objects, but exchange mass with the medium in 
which they are embedded. This poses a general prob- 
lem and uncertainty in defining the SFE. Imagine, for 
instance, a cloud accreting more gas from the outside 
than is turned into stars inside in the same amount of 
time. In that case, the SFE would go down, even though 
stars form continuously, simply because the total gas con- 
tent in the region of interest increases due to accretion 
from the outside. An even more remarkable thought ex- 
periment is one in which star formation has stopped at 
a given SFE in a given volume, and suppose the stars 
formed in that region now blow away all the remaining 
gas in their surrounding and drive it out of that volume 
(e.g., by driving winds). According to the standard defi- 
nition of the SFE in Equation would remain the 
same in this process, but Mgas would go down, leading 
to an increasing SFE, even though no more stars form. 
This is an example emphasizing that one must be quite 
careful with interpreting SFEs in open systems. Such 
definitions of the SFE in open syst ems are discussed in 
Vazquez-Semadeni et al. (2010) and Feldmann & Gnedin 



Shirley et al.|[2QQQl IKainulainen et"!!! [20091 [AF 



zoumanian et al.||2011| [Schneider et al.||2012| : 



( [2011[ ). 



9. SUMMARY AND CONCLUSIONS 



We studied the star formation efficiency (SFE) in simu- 
lations of supersonic, self-gravitating, magnetized turbu- 
lence. In Paper I, we found that the simulations yield star 
formation rates (SFRs) varying by orders of magnitude, 
depending on the virial parameter, on the forcing of the 
turbulence, and on the sonic and Alfven Mach numbers, 
consistent with observed SFR column densities. Here, 
we focused on the SFE and find that the density and 
column density distributions change significantly when 
star formation occurs and when the SFE increases. We 
find that the column density scaling is particularly sen- 
sitive to the SFE, enabling us to estimate the SFE from 
column density observations, without requiring a priori 
information on star formation activity or YSO counts. 
More detailed results are listed in the following. 

1. We find that star formation affects the column den- 
sity structure of clouds (see Figures [T| and [2). The 
effect of star formation on the global structure is 
stronger in low Mach number clouds and for clouds 
with solenoidal forcing than for clouds with high 
Mach number and/or compressive forcing, because 
strong turbulence alone produces high-density fila- 
ments, without requiring gravitational contraction. 

2. The volumetric and column density PDFs become 
wider with increasing Mach number and more com- 
pressive forcing (see Figures [s] and [4|. In con- 
trast, increasing the magnetic neld strength yields 
smaller standard deviations of the PDF (see Fig- 
ure [5|, as expected from Equation ([5|. 

3. The PDFs of our simulations with extremely dif- 
ferent forcing, sonic Mach number, and magnetic 
field strengths, all develop power-law tails of flat- 
tening slope with increasing SFE. We find that the 
high-density tails of the PDFs are consistent with 
equivalent radial density profiles, p oc r~'^ with 
K = 1.5-2.5 depending on SF E, consistent with 
column density observations (Motte et al. 1998 



4. The velocity spectra from our simulations with dif- 
ferent forcing, sonic Mach number, and magnetic 
field are all consistent with the scaling inferred 
in molecular cloud observations, v oc 

^1/2 (Larson! 

19811 IS olomon et al.[[1987[ [Falgarone et al. 19W 
Hssenko pf Mac Low||2Q02| [Heyer Brunt||2004 
Roman-Duval et al.[[2011), suggesting highly com 



pressible, supersonic turbulence of Burgers type on 
large scales. We find evidence for a transition to 
subsonic, Kolmogorov turbulence with v cx ^^/^ on 
scales smaller than the sonic scale (see Figure [6|. 
The velocity spectra are largely insensitive to star 
formation. 

5. Unlike the velocity spectra, we find a strong change 
of the density and column density spectra with the 
onset of star formation. The slope a of P3d(p, ^) 
/c" or equivalently of the A-variance spectrum of 
column density, cr^(I],^) cx which are both 
related to the column density spectrum P2d(S,^) 
via Equations ([isj and ( [2Q| , switches sign from 
a < for SFE to a >~0 when star formation 
proceeds, SFE > (see Figures [7|and[8]). 

6. The change in slope a occurs in all our models, 
despite the wide range of physical cloud condi- 
tions probed in the simulations {M ~ 3-50 with 
solenoidal, mixed, and compressive forcing of the 
turbulence, and magnetic fields ranging from the 
super-Alfvenic to the trans- Alfvenic regime; see 
Table [T]). We confirmed this result for the slope 
a with an independent measurement based on the 
A-variance technique, suitable for molecular cloud 
observations (see Figures [9]and[l0|. 



7. Studying the detailed dependence of the slope a on 
SFE in all our numerical models, we find that the 
SFE can be estimated by measuring a (see Fig- 



ure 



11). We provide a fitted relation, SFE (a) in 



Equation (22) for each numerical model set with 



best-fit parameters listed in Table [2] 

8. We compared observational measurements of a in 
an H I cloud, in several giant molecular clouds and 
sub-regions, and in a dense star- forming core, with 
values ranging from a = —1.6 to +1.6, all hsted 
in Table [Sj We infer the SFE for each cloud or 
region studied and show that negative a typically 
indicates SFE ~ 0, consistent with observations. 
Clouds with a ~ and larger do form stars and 
exhibit a range of SFEs, which is monotonically 
increasing with a. Estimates of SFE from a are 
consistent with independent measurements based 
on YSO counts, where available. However, the 
SFE(a) relations found here provide an indepen- 
dent estimate of SFE based on a column density 
map alone, without requiring a priori knowledge of 
star-formation activity or YSO counts for a given 
cloud. 

The overall agreement between our simulations and ob- 
servations of Galactic clouds concerning the SFE studied 
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here and the SFR in Paper I is promising. We conclude 
that supersonic, magnetized turbulence is a key process, 
likely controlling the SFRs and SFEs of molecular clouds 
in the Milky Way and potentially in other galaxies. 
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APPENDIX 

A. DEPENDENCE OF THE DENSITY PDFS 
ON RESOLUTION 



B. DEPENDENCE OF THE FOURIER 
SPECTRA ON RESOLUTION 



Figure [12] shows a resolution study of the density PDF 
p{s) for our two Mach 10 models with solenoidal and 
compressive forcing (left and middle panel), and for our 
Mach 50 model with compressive forcing (right panel) 
at SFE = and SFE = 20%. We compare resolutions 
A^res = 128, 256, and 512 for the Mach 10 models. We see 
no significant resolution dependence for A'res ^ 256 in the 
log- normal part of the PDF, but a clear resolution depen- 
dence in the power-law tail, stretching the PDF to higher 
s when the resolution is increased. In the solenoidal forc- 
ing case (left panel), we find a break in the power law at 
Sbreak ^ 4.0, 5.2, and 6.5 for A^res = 128, 256, and 512, 
respectively. These breaks seem to correlate with the 
density threshold for sink particle creation. Equation ([T]), 
^sink = ln(psink/po) ~ 3.1, 4.5, and 5.9, shown as vertical 
short lines on top of the abscissa {N^es = 128: dotted; 
256: dashed, and 512: solid), but are slightly higher than 
that. The PDF for densities s < Sbreak seems converged, 
while the flatter part for which s > Sbreak shifts to higher 
densities when the resolution is increased. The compres- 
sive forcing model with Mach 10 (middle panel) does not 
show such a break, but a cutoff at large 5, which shifts to 
higher s with increasing resolution, similar to the cutoff 
seen in the solenoidal forcing model at large s. The PDF 
of the Mach 50, compressive forcing case (right panel) 
does not show a significant resolution dependence for the 
range of densities and Ps analyzed here {N^es = 256: dot- 
ted; 512: dashed; and 1024: solid), even though the sink 
particle creation densities are quite low (see the vertical 
lines on top of the abscissa). 



Figure [13] shows the same resolution study as in Fig- 
ure [l2j but for the density and column density power 
spectra as in Figure [?[ This resolution study demon- 
strates that density and column density spectra are con- 
verged for k < A^res/15. In contrast, velocity spectra are 
only c onverged for k < A^res/30 (Federrath et al.| |2010b[ 
2011c). This is because resolving the kinetic energy con- 
tent of a vortex (which is a vector quantity) requires at 
least 30 grid cells. Here we find that about 15 grid cells 
are sufficient to resolve turbulent and collapsing density 
structures. This in turn suggests that resolving the Jeans 
length with 15 grid cells would be sufficient to capture the 
density structures in collapse simulations, but is insuffi- 
cient to resolve kinetic vorticity and potential small-scale 
dynamo amplification of magnetic fields (both of which 
are vector quantities), which still require s at least 30 grid 
cells p er Jeans length ('Sur et al.||2010l Federrath et al. 
2011c). The slopes a that we estimate in Figures [7}|TU 
are obtained in a converged range of scales, and hence 
these slopes are not affected by numerical resolution. 



C. 



DEPENDENCE OF THE DENSITY PDFS 
AND FOURIER SPECTRA ON THE 
RANDOM SEED 



Figure 14 (left) indicates statistical fluctuations mostly 
in the low- and high-density tails of the PDFs when 
different random seeds of the turbulence (seedl, seed2, 
seed3) are considered. This gives us an impression of 
the statistical uncertainties of the PDFs. The right 
panel shows the density and column density spectra of 
the same models, also indicating statistical fluctuations, 
which manifest in different amplitudes of the spectra, 
but the power-law slopes do not vary significantly with 
random seed. 
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